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RESEARCH SUMMARY 


This report describes a user-'Oriented computer system that 
allows fire managers to quickly and easily analyze climatological 
data for the purpose of predicting the probable occurrence of 
desired conditions for prescribed fire. Two separate computer pro- 
grams, RXWTHR and RXBURN, make up the system. Both programs are 
designed to use National Fire Weather Data Library climatological 
data. Among the features of these programs are: adjustment of 
fuel moistures from the weather station site to the fire site; a 
newly developed duff moisture model; and the capability to simul- 
taneously consider up to 15 prescription factors. 

Step-by-step instructions for use of the programs are given 
and demonstration runs are provided. Documentation of all pro- 
grams is included. A computer terminal operator's manual for 
programs RXWTHR and RXBURN has been published as a separate Gen- 
eral Technical Report (Bradshaw and Fischer 1981) , 
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INTRODUCTION 


The timely occurrence of desired burning conditions has been and continues to be a major 
source of uncertainty associated with using fire as a resource management tool. Although such 
uncertainty cannot be eliminated, it can be reduced by using information contained in clima- 
tological records (Fischer 1978); for example, analysis of climatological data can describe an 
areals average weather and also likely deviations from this average situation. A fire manager 
can use such information to identify periods when desired burning conditions might occur, 

Further analysis can result in calculation of the probability of desired fire conditions 
occurring during the periods identified (Furman 1979). 

The value of climatological analysis as an aid for scheduling prescribed fires can be 
summarized as follows: 

1. Prescribed *fire schedules can be geared to the number of probable burning days. Over- 
scheduling of fires, especially during traditional burning periods, can be minimized. 

2. Prescribed fire opportunities outside traditional burning periods can be identified 
insofar as available weather data will allow. 

3. Project priorities can be more nearly matched to the probable occurrence of prescribed 
conditions. Some fires require a narrow range of burning conditions for successful accomplish- 
ment. These fires can be given top priority for scheduling during periods when such conditions 
are most likely to occur. 

4. Fire management prescriptions for preplanned use of chance wildfires can be evaluated 
before they are implemented. The frequency and time of occurrence of prescribed fire conditions 
can be identified. This information can then be used to evaluate the probable impact of the 
prescription and then compared with the probable performance of alternative prescriptions. 

Despite its obvious value, fire managers rarely have had the means to easily analyze 
climatological data when planning fire use. This is due, in part at least, to the tedious 
nature of the task. It is also due, in part, to weather data problems: missing data, lack 
of data during spring and fall, and, in some cases, unreliable data resulting from poor 
observation and poor weather station maintenance. 

This report describes a user-oriented computer system that allows fire managers to quickly 
and easily analyze climatological data for the purpose of predicting the probable occurrence of 
desired prescribed fire conditions. Instructions are provided to assist both the fire manager 
and the computer specialist through the various steps required to use the system. Finally, 
documentation of the computer programs used is provided, A computer terminal operator's manual 
has been published as a separate report (Bradshaw and Fischer 1981). 

SYSTEM DESCRIPTION 


Two separate, but interrelated computer programs, RXWTIIR (Prescribed Fire Weather) and 
RXBURN (Prescribed Fire Conditions) make up the system. Both programs are designed to use the 
climatological data base stored in the National Fire Weather Data Library (NFWDL) (Furman and 
Brink 1975), which is located at the USDA Forest Service's Fort Collins Computer Center (FCCC) . 
RXWTIIR and RXBURN are stored in the Region 1 program library at FCCC and are available via remote 
terminal or batch processing to all who have access to FCCC facilities. 

Information Provided by the programs 

Program RXWTHR provides climatological summaries and co-occurrence frequencies of user- 
selected fire weather and fire danger rating parameters. The 16 parameters available for 
analysis are: 
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1. State of the weather 

2 . Temperature 

3. Relative humidity 

4. Wind direction 

5. Wind speed 

6. Maximum temperature (last 24 hours) 

7* Minimum temperature (last 24 hours) 

8. Maximum relative humidity (last 24 hours) 

9, Minimum relative humidity (last 24 hours) 

10. Precipitation duration 

11. Precipitation amount (last 24 hours) 

12. 1 hour fuel moisture 

13. 10 hour fuel moisture 

14. 1978 National fire Danger Rating System (NFDRS) Energy Release Component (ERC) 
(with selected fuel models) 

15. NFDRS Burning Index (BI) (with selected fuel models) 

16. Duff Moisture (24-hour average for entire layer) 

The user is reminded that Program RXWTHR can be used to summarize and analyze only those 
parameters that are included in NFWDL for the weather station of interest. The importance of 
maintaining complete records at the NFWDL and keeping them current tlu'ough the Adminstrative 
and Forest Fire Information Retrieval and Management System (AFFIRMS) (Helfman and others 1975) 
cannot be overly stressed. 

Examples of the type of information provided by RXTHWR are shown in figures 1^ 2, 3, and 4, 
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Figure 1 . --Example of Program RXWTHR summary table for temperature. 
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Figure 2. --Example of Program RXWTHR summary table for wind direction. 
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Figure 4. --Example of Program RXWTllR three-way co-occurrence table for temperature^ relative 

humidity, and windspeed. 

Program RXBURN provides analyses of local prescription condition frequencies based on 
user-defined prescription conditions* The 16 parameters available for use as prescription 
conditions are: 

1. State of the weather 

2. Temperature 

3. Relative humidity 

4. Wind direction 

5. Wind speed 
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0. MtiMinuTii temperature (last 24 hours J 
7. Minimum temperature (last 24 hours) 

S» MaAimum relative humidity (last 24 hours) 

9. Minimum relative humidity (last 24 hours) 

10. ^Number of days since last measurable precipitation 

11. Precipitation amount (last 24 hours) 

12. 1 hour fuel moisture 

15. 10 hour fuel moisture 

14. 1978 NFDliS ERC (with selected NFDRS fuel models) 

15. 1978 NFDRS BI (with selected NFDRS fuel models) 

16. Duff moisture (24-hour average for entire layer). 

As many as 15 of these parameters may be used in a single prescription. Information pro- 
vided by RXBURN is exemplified in figures 5, 6j 7, and 8, 
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(PERCENT) 

month of highest prescrition frequency 
(Percent probability) 

10 DAY PERIOD O'' HIGHEST RX FREQUENCY BEGINS 
(PERCENT PROBABILITY) 


PREFERABLE 

ACCEPTABLE 

UNACCEPTABLE 

a 

51 

SO 

Sfc 

37* 

5B* 

SEP 

JUL 

OCT 

RIv 

45* 

71* 

JJN 1 

OCT 1 

JUL 1 

MAY 1 

\f)% 

57* 

7fl* 


Figure S.— Example of Program RXBURN summary table. 
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PROGRAM RXBU^NI RUN NO. 1 DISTRICT! FIRELAB 


FOREST! MONTANA 


PAGE NO. ^ 


Prescription occurrence ry lo day period and month 




tt 

* PERIOD *» PREFERABLE OAVS *• ACCEPTABLE DATS »• UNACCEPTABLE DAYS 

• MONTH BEGINS NO. DAYS MEAN NUMBER PERCENT MEAN NUMBER PERCENT MEAN NUMBER PERCENT 


tf 


ft 

HAY 

1 

104 

0 

5 

S% 

ft 

HAY 

11 

109 

0 

5 

S% 

ft 

HAY 

21 

118 

1 

B 

73i 

ft 

MAY 

total 

331 

1 

IB 

59i 


« 


1 

IB 

17% 

5 

81 

7B% 

2 

33 

30% 

4 

71 

SB% 

3 

AO 

34% 

4 

70 

59% 

6 

91 

27% 

14 

222 

fe7% 




ft 

JUN 

1 

llO 

1 

11 

10% 

3 

50 

A5* 

3 

49 

45% 

ft 

ft 

JUN 

U 

109 

0 

4 

4% 

2 

39 

36% 

4 

66 

61% 

ft- 

ft 

JUN 

21 

10ft 

0 

6 

6% 

3 

45 

A2« 

3 

55 

52% 

ft 

ft 

JUN 

total 

325 

1 

21 

6% 

8 

134 

Al« 

11 

170 

52% 

ft 


« » 


ft 

JUU 

1 

i3a 

0 

3 

2% 

5 

78 

S7« 

4 

57 

41% 

ft 

ft 

JUL 

11 

150 

1 

6 

5% 

4 

62 

AH 

5 

80 

53% 

ft 

ft 

JUL 

21 

165 

0 

5 

3% 

4 

63 

38« 

6 

97 

59% 

ft 

ft 

JUL 

total 

453 

1 

16 

4% 

13 

203 

A5» 

15 

234 

52% 

ft 


« tt 


»«#«»!» ft im «»»#«»*«*«#»««*««»»»»»«•«»«»«««««•«»«»«»»« «#»»»»•»» 


ft 

AUG 

1 

159 

0 

7 

4% 

4 

66 


5 

86 

54% 

ft 

ft 

AUG 

n 

146 

0 

3 

2% 

4 

56 

38« 

5 

87 

60% 

II 

ft 

AU3 

21 

16o 

1 

10 

6% 

4 

67 

A2* 

5 

83 

52% 

ft 

ft 

AUG 

TOTAL 

465 

1 

20 

4% 

12 

IB9 

AH 

16 

256 

55% 

ft 


ft ft 

ftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftfttfftftftftftftftftftftftftftftftftftftftftftftftftftftftttftftftftftftftttftftftftftftftftftttfttftftftftftftft 


ft 

SE» 

1 

137 

1 

12 

9% 

3 

55 

40% 

4 

70 

51% ** 

ft 

SEP 

11 

129 

1 

11 

9% 

3 

45 

35% 

5 

73 

57% » 

ft 

SEP 

21 

117 

1 

11 

9% 

2 

39 

33% 

4 

67 

57% * 

ft 

AA 

SEP 

TOTAL 

383 

2 

34 

9% 

9 

139 

36% 

13 

210 

5S% » 

ft ftftftftftftftft ftft ft « ftft ft ft ftft ft ft ttftfttft ft ftftftttftftftftftllftftft ftftftftftftftft ffttftftftftftftftftftftftft ftftftftft 

ftft ftft ftft ftftftftftftftftftftftftftftftftft ft 

ft 

OCT 

1 

90 

1 

9 

10% 

1 

20 

22% 

4 

61 

68% * 

ft 

OCT 

11 

90 

0 

3 

3% 

1 

21 

23% 

4 

66 

73% » 

ft 

OCT 

21 

70 

0 

0 

0% 

1 

22 

2B% 

4 

56 

72% » 

ft 

ft 

OCT 

total 

258 

1 

12 

5% 

4 

63 

24% 

11 

183 

71% » 


ft lift ftftftftftftftft ft lift ft ftftftftftftftftftftHftft ftftftftftftftft ft ftft ftftftftftttftftftftft ftft ftftftftft ft ft 

ft 

ft 

NOV 

total 

0 

0 

0 

0% 

0 

0 

0% 

0 

0 

0% • 

ft 

«l 

ft 

ift ftftftft fti 

TOTAL 

IftftftftftftftI 

TOTAL 

221s 

Iftftftftft 

B 

ftftftftftftft 

121 

ftftftftftftftfti 

5% 

► ftftftftft! 

51 

iftftftftftfti 

819 

37« 

80 

1275 

58% * 


ft ft 


ft ftft ft<ft»ft»ftftftft ftftftftft ftft ftft ftftftftftftftft ftftftftftftftftftftftftftftftftftftftftfttftftttftftftftftftftftftftftftftftftftftftftftftftftftftftftftft ftft ftftftftft 


Figure 6. — Example of program RXBURN prescription frequency of occurrence table^ 
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PPOS'JtH «XBU>NI ’UN NO. 


1 DISTRICTI FIRELAR 


'■OREST I HOAdANA 


PAGE NO, ) 


PRESCRIPTION RUN UFNGTH SUMMARY 




PERroO 


•• PREPEffABlE DAY RUNS 
* PERCENTILES ♦ 



MONTH 

REOrNS 





HAlf 

1 


MAT 

11 


M4r 

21 


HI4r 

TOTAL 





JUN 

1 


UJN 

11 


JJV 

21 


JJN 

TOTAL 





JJ. 

1 


JUw 

n 


JUl 

21 


UJl 

TOTAL 


** ACCEPTABLE DAY RUNS »» UNACCEPTABLE OAY RUNS » 
PERCENTILES ^ » PERCENTILES » » 

I MEAN ?5TN MEDIAN 7STH MEAN ?5TH MEDIAN 7HTH ^ 



Figur. 7,-.Exa.ple of Pr„gra„ „„„ |n-oscri|.. i„„ r„„ lc„|.th 
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program rxbu^n: 


1 


DISTRICT: FIRElAR 


^ORESTS M3NTAtgAi 


PAGE NO* A 



PR0B^3ILI rv 

Of MFETING 

PRESCRIPTION 2* A\»D 

3 DAYS IN THE FUTURE 




0 0000 

« 





MONTH: 

MAY 





0 

tt 

todays 

00 0 

TOMORROW 

^*00 

0000 

? DAYS 00*0 

0000 

3 DAYS 

000 0 

M DAYS 

0 


cond 

PREF 

ACCP 

UNAC 

PREE 

ACCP UNAC 

PRFF 

ACCP 

JNAC 


0 

0 

0 

^REF 

17% 

50% 

33% 

9% 

3B% S3% 

7 % 

32% 

61% 

P9EF 18 

0 

0 

0 

ACCP 

12% 


AU 

8% 

35% S6% 

7 % 

31% 

62% 

ACCP 85 

0 

0 

0 

JNAC 

2% 

18% 

80% 

A% 

2A% 72% 

5% 

26% 

6B% 

UMAC 218 

0 

n 

000000000«0000000000000»00»0«000000000000000tf0000000tfa«000«00000000«000000-0000000000000000 



0 





month: 

JUN 





0 

0 

todays 

00 0 

tomorrow 

000 

0000 

2 days 

0000 

3 DAYS 

000# 

M DAYS 

0 

0 

COMO 

PREF 

ACCP 

UNAC 

PREF 

ACCP UNAC 

prFF 

ACCP 

JNAC 


0 

0 

0 

3REF 

20% 

60% 

20% 

10% 

5H 39% 

B% 

45% 

47% 

PRFF 20 

0 

0 

ACCP 

9% 

56% 

35% 

8% 

A7% A5% 

7% 

4 3% 

49% 

ACCP 132 

0 

0 

0 

JNAC 

3% 

PB% 

69% 

5% 

37% 58% 

6% 

40% 

54* 

UMAC 173 

0 

H 



0000000000000000 

...» 

0 





MONTH? 

JUl 





0 

n 

0 

todays 

»»» T0^^0RR0W 

0 00 

000 0 

2 DAYS 

000 0 

3 DAYS 

0 00 0 

N DAYS 

0 

0 

cond 

PREF 

ACCP 

UNAC 

PREF 

ACCP UNAC 

PREF 

ACCP 

JNAC 


0 

n 

« 

^REF 

0% 

75% 

25% 

A% 

48% 48% 

4% 

45% 

51* 

PREF 16 

0 

0 

0 

0 

ACCP 


52% 


A% 

46% 51% 

4% 

45% 

52* 

ACCP 206 

0 

1^ 

0 

JNAC 

3% 

36% 

61% 

3% 

A3% 54% 

4% 

44% 

52* 

UNAC 227 

* 


# 


Figure 8. --Example of program RXBUUN prescription occurrence contingency table. 


Information Needed to Use the Programs 

Both programs require climatological data from the National Fire Weather Data Library, 
Information on the NFWDL data base structure and use is detailed by Furman and Brink (1975), 

For completeness, instructions for accessing NFWDL data for use in RXWTHR and RXBURN are outlined 
in the User's Guide in the next section of this report. 

The user must also furnish information requested on the information worksheets For RXWTIIR 
and RXBURN (figs. 9 and 10). 
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P ro — 


^;re .vojtiier at Lon [ntoniJt ion Station name_ 


iz I evati^n 


tt .Latitude 


, Climate class— 


, Slope c lass— ^ 


Ffiel model- 
'i ear be^i n 


.Last trost ^ 
,^ear end 


.Grass type : annual 


.Date begin^ 


. perennial 


.Date end 


Site Ad’Ubtment Factors tif anvl 


Aspect f l=north, H-east, S^south, d^west ). Site elevation^ 

Canopv cover I l=open, 2-closed ] 

Duff/Soil Horizon Intormatian if Duff Moisture selected. See User's Guide, appendix D) : 
La>er Duff^ioil Type Thickness ♦*»*iri*»********#*^Ht*^*******##**i»*******^****^*,»*iri*i 


rt 


I 

cm 

* 

* 

y 

Use up to 80 characters 


cm 

■* 

w 

w 

:/ 

See User’s Guide, append 

3 

cm 

★ 

■X 

1/ 

See User’s Guide, append 

■4 

cm 

■X 

♦ 

£/ 

For NFDRS indices only. 



» 


append i.x A 

S 

cm 

* 

5/ 

For NFDRS indices only 


Summary Tablet sj Requested (Select up to 15 ] : 

State of the v^eather 

^Temperature (degrees F) 

Relative humidity {"!>) 

^'^ind direction {S point) 


Jund speed tmi/h) 

temperature Cd h, deg. Fj 
_Min cemnerature \Z^ h, deg. FI 
_Ma\ relative humiditv [2^ h, M 


^^lin relative humidity (24 h. '^0 

Precip duration (last 24 h) 

Precip amount (24 h, 0.01 in) 

1 hour fuel moisture (‘o) 

10 hour fuel moisture {%) 

^MFDRS ERC 

^NTDRS BI 

Duff Moisture (24 h avervige, 


Tablelsl Req uested ilf 2-hav table desired leave last 


selec ted, Wind Direction must always be listed fir st.-): 
_vith 
with 


space blank. If 


with 


with 


-I 


jvith_ 

_with_ 

with 


with 


■ with 
with 


Figure 9.— RXimiR user information sheet. 
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-- UOCK lUrN -^ntic 1 


Total Mumber Ditferent Stations in This Run 

User's Name .Subunit .iJnit 

Pro j 


Fire Weather Station Information : Station name__ ,No. 

Elevation^ ft . > Latitude .Climate cLass -"^ .Slope class^^^ 

Fuel modeli^^ .Last frost^;|||^ .Grass type^'^: annual .perennial 

\ear begin^j .'tear end .Date begin 

Site Adiustment Factors (if any) - 

Aspect { l=north, 2=east, 5-south, 4=west ), Site elevation ■ 

Canopy cover ( l^open, 2=closed ) 


Duff/Soil Horicon Information (if Duff Moisture selected. See User^s Guide, \ppendix D) 
Layer Duff/Soll Type Thickness 


★ 


1 

cm 

* 

* 

y 

Use up to 80 characters 

') 

cm 

★ 

k 

k 

y 

See User's Guide, appendix C 

5 

cm 

k 

k 

y 

See User's Guide, appendi.x B 

4 

cm 

k 

k 

K 

4/ 

For NFOP.S indicies only. See 

5 

cm 

k 

k 

k 

y 

appendix A 

For .vFDRS indices only 


See User ' s Guide , 


Prescription Factor Selections (Check and set limits for up to 15 factors. 1 : 

Preferable Rx Acceptable Rx 

Factor Minimum Maximum Minimum Maximv 

^State of the weather 

^Temperature (deg. F) 

^Relative humidity (“■) 

Wind direction (S point) 

Wind speed (mi/h) 

^^lax temperature (24 h, deg. F) 

^Slin temperature (24 h, deg. F) 

Nlax relative humidity (24 h, %) 

>tin relative humidity (24 h, "o) 

^Days since last precipitation 

Precip amount (24 h, 0.01 in) 

hour fuel moisture (%) 

10 hour fuel moisture (%) 

^NFDRS ERG 

NFDRS 

Duff moisture (24 h average, "i) 


Figure 10.— RXBURN user information sheet. 


Program Use 

Ideally, the programs should be used sequentially. RXWTHR will give the manager a method 
for screening possible combinations of prescription factors to eliminate those which have a 
low probabilitv of occurring simultaneously, as well as supplying a local fire climatological 

MHiiiiury, 

Once a feeling for the occurrence and timing of desired prescription conditions is reached, 
gXBURN may be used to provide a detailed summary of prescription occurrence frequencies, 

Because of the costs involved, RXWTHR should be run as few times as possible for any given 
station, while RXBURN may be run many times for different prescriptions and planning processes. 
Hach program is, however, completely independent of the other and sequential use is not 
mandatory. 


Options and Features of the Programs 

f'ORTRAN IV code is ASCII standard, 

1978 National Fire Danger Rating System (NFDRS) algorithms and subroutines are used by 
the programs where applicable (Deeming and others 1977). 

The routines operate on simple input streams, 

The programs contain input error checks that print error messages. 

The programs offer 1978 NFDRS indices, Energy Release Component (ERG) and Burning Index 
(BI) as prescription parameters based on NFDRS fuel models. 

The programs contain site adjustment factors to adjust fuel moisture values (and 
resulting ERC^s and BI's) to locations with elevations, aspects, and canopy covers different 
than those of the user-specified station, 

RXbrilR allows as many as 15 summary tables and five co-occurrence tables to be produced 
in a single run. 

^ RXDURN allows as many as 15 prescription factors to be simultaneously considered in a 
prescription condition. 


As many as 99 fire weather stations may be analyzed for summaries or proscription 
occurrence in a single run. ^ ^ 


Repeated runs 
ditions to be analyzed 
first rim computations. 


of the same fire weather station allowing different prescription con- 
in one run is possible through the construction of a disk file to hold 


A newly developed duff moisture model is available (Fosberg 1975). 


factor, 

burning 


RXBURN allows two prescription condition 
One is an ideal or ’’preferable'* burning 
condition. 


ranges to be entered for each prescription 
condition, the other is an "acceptable” 


RXBURN allows the number of days since last precipitation as 


a prescription condition. 
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System Limitations 

Programs RXWTHR and RXBURN are designed to use archival data in the National Pire Weather 
Data Library and Brink 1975). Gibson [1977) and others l\ave noted a limitation in the 

use of National Fire Weather Data Library records for planning prescribed fire- -observations are 
taken only once a day over the length of the fire season. This results in two constraints on any 
system that uses the data. First > the single observation must be considered to be representative 
of the entire day, when in fact it represents the ’’typical worst" conditions [midafternoon, 
southwest aspect, inidslope, open canopy). Second, weather observations are rarely recorded for 
more than 5 months each year during the "fire season," Data therefore tend to be scarce in 
prefire and postfire season periods- -times that prescribed fire operations are generally most 
active. 

To reduce computer storage requirements, program RXWTHR analyzes a maximum of 5 months of 
weather data from the National Fire Weather Data Library in a single run. This is usually not a 
limitation of the system since few stations have more than 5 months of data on file for each 
year, For those stations with more than 5 months of fire weather data, two or three runs of 
RXWTHR may be needed to obtain the desired climatological summaries, The 5-montli restriction 
does not apply to program RXBURN. 

A newly developed duff moisture model is available for test use (Fosberg 1975), The model 
IS based on theoretical considerations of water transport through the litter, duff, and soil 
horizons as affected by daily weather observations. The model is still in the validation stage, 
and results should be viewed as tentative and used with caution. Preliminary results indicate 
the model will consistently underestimate duff moisture content. 

USER’S GUIDE 


Step-by-step instructions for using RXWTHR and RXBURN follow. Demonstration runs of both 
programs are included to support these instructions. The instructions are written so that anyone 
familiar with the use of remote tei'minal or batch processing equipment should be able to success** 
fully run the programs. 

In organizations that employ computer technicians, the fire manager’s role In running these 
programs may be limited to completing the user information sheets (figs. 9 and 10). 

Create a Data File 

Two items of information are needed to obtain data from the NFWDL. First is the six-digit 
code [or codes) of the fire weather station(s) to be analyzed and the years of data to be ana- 
lyzed. Second is the file name in the library that contains the lowest station code that will be 
used in the analysis. For example, if the stations to be analyzed are 034567, 245789, and 
003452, only the file name that contains station 003452 is needed. 

Obtain a File Name 

There is a possibility that a local fire weather officer or computer specialist has a 
current listing of the files in the library. If not, obtain a file list in the following manner: 

1 2 3 4 5 

12345678901234567890123456789012345678901234567890 


@RUN, . . . 

@aSG,A PIREDATALIB*PROGRAMS. 

OXQT FIREDATALIB*PR0GRAMS.L1STFILES 
@FIN 
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System software will then respond with a listing of file names and stations in the file. An 
evampie of the general format follows. (Note: ssssss represents the six-digit station code, and 
vy represents the last two digits of the year that data begins (FROM) or ends (THROUGH); nn, mm, 
and 00 represent assorted numbers and letters of file names.) 


FILE UmE 

STATION 

YEAR LIMITS 

DATE OF LAST 


FROM 

through 

UPDATE 

riREDATALIB*nn“inm 

ssssssyy 

ssssssyy 

mmddyy 

FIREDATALIB*00 

ssssssyy 

ssssssyy 

mmddyy 


Scan the station year limits column until the group containing the lowest six-digit station 
code of stations to be analyzed is found. The entire file name (FIREDATALIB*mm-nn) is to be used 
in place of ’’FILE'* in the following data acquisition sequence. 

Create a Card Image File 

In creating a card image file for use in RXWTHR or RXBURN, it is wise to generate a user 
program file of the data. This allows the data file to be stored in the Mass File Directory at 
rCCC for 6 days from the date of creation, allowing subsequent runs of RXWTIIR or RXBURN to access 
the data with recreation of the file on each run. This is particularly helpful in the event of 
input stream errors resulting in program termination prior to complete execution. 

This process is accomplished by executing the following control sequence at the FCCC, 

(Note: ^galn, ssssss represents the six-digit station code; yy the year to begin data Inclusion 
in the file; and the second yy the year to end data inclusion. If all available years are 
requested, use yy = 00 and yy - 99.) 

12 3 4 

1234567890123456789012345678901234567890 

ORUN, ... 

QASG,A F I REDATAL IB* PROGRAMS. 

QaSG,a file. 

QUSE 2., FILE. 

0ASG,UP NAME. (NAME MAY BE ANY NAME) 

0USE 15. , NAI’lE. (USE SAME NAME AS PREVIOUS LINE) 

^XQT FIREDATALIB*PR0GRAMS.GETDATA2 


ssssssyy ssssssyy 

etc . f until all stations are listed in ascending order. 
@EOF 


OFREE 2, 

aASG,T 2, ASSIGNS TEMPORARY DISK FILE USED IN PROGRAMS) 


1 ages in each station, Tlie file is now ready for analysis by nrogram RXimn 
-ead weatlier data from logical unit 15. ^ ^ program RXWTHF 


the number of card 
or RXBURN, wliich 
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Subsequent runs using the same station within the next 6 days would use the following 
sequence in preparing the data file for program use: 

1 2 
12345678901234567890 

Orun, . . . 

@ASi3,A NAME. (NAME IS ONE PROM ORIGINAL RUN) 

Ouse 15. ,name. 

@ASG,T 2. 


Using Program RXWTHR 

To use RXWTHR, four steps must be followed, 

Step 1, Pill out user information sheet (s), (See fig, 9,) 

One sheet should be completed for each station to be summarized, or for each summary for 
a station. The entries on the sheet are mostly self-explanatory, and are briefly summarized 
here for completeness. 

a. Total number of stations in run . --This value is the number of different stations that 
will have summaries produced in one run of RXWTHR, 

b. Fuel model, grass type, and date of last spring frost . --These values need only be 
entered when NFDRS indices are requested in output tables, NFDRS fuel model descriptions are 
detailed in appendix A. 

Slope class , --This value must be entered if NFDRS indices or site adjustment options 
are selected. Slope classes are defined in appendix B, 

d, NFDRS climate class , --This value should be entered on all information sheets and as 
detailed in appendix C, 

0 , Activity information . --User comments and documentation are entered here. As many as 
80 characters are allowed, Suggested entries include user name, project objectives, site 
adjustment factors, district, and other information. 

f, Duff/soil horizon information . --This information is only required by the selection 
of duff moisture as an output parameter. Instructions are detailed in appendix D. 

g, Summary tables . --If summary tables are desired, select parameters for inclusion. 

h, Co-occurronce tables . --If co-occurrence tables are desired, select up to five by 
filling the blanks with the desired parameters from the total list in the preceding summary 
section. The only limitation is that wind direction, if selected, must be the first param- 
eter entered in a table composition. 

Step 2, Create data file for program use and load program RXWTHR, 

File creation was previously covered; program loading is done by the execution of the 
following control cards: 

12 3 

123456789012345678901234567890 

@ASG,A CSSG*R1LIB. 

@XQT CSSG*R1LIB, RXWTHR 
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step 3* Transfer information from user sheet (s) to machine readable formats. 

The information from each information sheet constitutes a directive block that tells the 
program what options are requested and supplies other needed information. There must be one 
directive block for each station in the analysis and for each multiple run of a single station. 
Information entered in one directive block does not need repeating in subsequent station analysis 
nor multiple runs if the information does not change, Details of options and information entry are 
detailed below with examples. 

a. Total number of stations card . --This card must be the first in each directive block. 

The value is entered in columns 1 and 2, right justified, no decimal. 


1234567890 

nn 


b. Station, years, and dates of analysis cards . --The first card in this sequence instructs 
the program that three cards of station, year, and dates information follow. It consists of 
tl\e word STATION beginning in the first column. The first trailing card contains complete 
station information as detailed in table 1, This is then followed by two more trailing cards 
with additional information. The first identifies the years to begin and end data inclusion 
to the analysis, the second identifies the seasonal dates of data inclusion. If all available 
years are desired, use 1900 to 1999. Program RXIVTHR is restricted to 5 months of data analysis 
in a single run. For stations with more than 5 months of archival data, see the following 
section on System Options. For an example, consider Philipsburg Ranger Station for the sea- 
sonal dates of June 1 to September 15 for 1960 to 1977. 


1 2 3 4 5 6 7 

1234567890123456789012345678901234567890123456789012345678901234567890 


STATION 

PHILIPSBURG RS 243002 5280 46 H 3 P 3 0615 


T 3 6200 1 


YEARS 

1960 

1977 

DATES 

0601 

0915 


c. Activity information card . --The first card in this sequence tells the program that 
one card of activity information follows. It consists of the word ACTIVITY beginning in 
column 1. It is followed by a card that contains the information under activity on the user 
information sheet. 


123456 7 8 

12345678901234567890123456789012345678901234567890123456789012345678901234567890 

ACTIVITY 

FIRE WEATHER SUMMARIES FOR PHILIPSBURG, LARRY JONES, PROJECT 5 
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Table I » --station information card 


Input variable 

Limits of 
value ranges 

Card columns 
(inclusive) 

Format 

Information 

source 

Station name 

15 characters 

1-15 

4A4 

user 

Station number (6 digits) 

0-999999 

17-22 

16 

user 

Station elevation (ft) 

0-99999 

24-28 

F5.0 

user 

Station latitude (degrees) 

0-99 

30-31 

12 

user 

NFDRS fuel model 

A-U (except M) 

33 

A1 

appendix A 

NFDRS slope class 

1-5 

35 

11 

appendix B 

Grass type 

A or P 

37 

A1 

appendix A 

NFDRS climate class 

1-4 

39 

11 

appendix C 

Date of last frost 

0101-1231 

41-44 

14 

user 

Repeat run ? (REPEAT)^ 

T or blank 

46 

LI 

user 

Save run ? (SAVE)^ 

T or blank 

48 

LI 

user 

Site adjustments ? (ADJUST)^ 

T or blank 

50 

LI 

user 

Site aspect code 

1-4 

52 

11 

user 

Site elevation (ft) 

0-99999 

54-58 

F5.0 

user 

Site canopy cover 

1-2 

60 

11 

user 

If the run is to save computations for the next run 
wise leave column 48 blank. 

to use, enter 

a T in card 

column 48, other 

If the run is to use computations 

from the last run 

, enter a T in 

card column 

46, otherwise 


leave column 46 blank. 

It the run is to make site adjustments, enter a T in card column 50. Be sure if site adjusnients 
are to be made that correct values are entered for aspect, elevation, canopy, and NFDRS slope 
class. 

If the run is not to make site adjustments, leave columns 50 through 60 blank. 


^Logical variables. 


d. Duff/soil information . --This input sequence begins with a card instructing the pro- 
gram as to how many layers (cards) of duff/soil information follow. It consists of the word 
DUFF in columns 1 through 4, and the total number of layers (2 to 5) in column 12, No decimal 
is punched. This card is then followed by as many as five cards that define the layer type 
and thickness (in centimeters) . Layer type is in columns 1 through 5 and thickness in columns 
6 through 10. Both values have decimals punched. 


1 2 
12345678901234567890 

DUFF 3 

1. 3.5 

2, 7.0 


7. 300.0 
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e. SuiTiniai^y table option cards , --This sequence begins with a card that instructs the 
pi'ogi^am that summary table options follow. The word SUMMARY is printed beginning in column 
1 and is followed by as many as 15 cards defining the parameters to be summarized. There is 
one card for each parameter checked on the sheet. The entire parameter name should be spelled 
correctly beginning in tlie sixth column of the card. Information in the parentlieses may be 
excluded. After the last option card has been printed, instruct the program that the last 
card lias been read by placing the word END beginning in the first column of the next card. 


12 3 

123456789012345678901234567890 

SUMMARY 

TEMPERATURE 

WIND SPEED 

END 


f. Co - occurrence table option cards . --This sequence begins by instructing the program 
that co-occurrence option cards follow. This is done by use of the word CO-OCCUR beginning 
in the first column. It is followed by as many as five option cards that define what tables are 
to be produced. Values are transferred from the user sheet to cards, one card per table. 

Make sure the spelling is the same as the parameters in the summary list, The first param- 
eter begins in the sixth column, the second in the 31st, and third (if any) begins in the 
56th column. Again, information in the parentheses should be excluded and the last option 
card is followed by an END card. 


1 2 3 4 5 6 7 

1234567890123456789012345678901234567890123456789012345678901234567890 


CO-OCCUR 

WIND DIRECTION 
TEMPERATURE 

END 


WIND speed 
RELATIVE HUMIDITY 


WIND SPEED 


g. The RUN card . --This card must be the last card in each directive block. It consists 
simply of the word RUN in columns 1 through 3 and signals that all cards in the block have been 
read. The program then begins execution. 


1 

1234567890 

RUN 


Step 4. Terminate program execution. 

The program will continue execution, reading one directive block at a time, and perform- 
ing the requested calculations until coming upon an end-of-file mark (@E0F) when reading the 
first card of the next directive block (Number of Stations). After all directive blocks have 
been entered, the program is terminated by the following commands; 
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1 

1234567890 

0EOF 

@FIN 


Program RXWTHR Error Messages 

Control sequence incorrect, program aborted. 

Cause: A control^ card was read that was uninterpretable by the program's internal diction- 
ary. Most common cause of an aborted program is misspelled control words or not 
starting the word in column 1 of the control card. Can also be caused by not 
placing an END card at the end of an option series. 

2. No summary or co-occurrence option specified, program aborted. 

Cause: There is neither a SUMMARY nor a CO-OCCUR control cai-d in the input stream. 

3. Station numbers on descriptor card and data file do not match, program aborted. 

Cause: Keypunch error on STATION card following STATION control card, or error in data 
file structure. Remember that stations must be analyzed in ascending order. 

4. At least one summary option card is incorrect, program aborted. 

Cause; One of the option cards trailing the SU^^MARY control card is misspelled or does 
not begin in column 6. 

5. At least one co-occurrence option card is incorrect, program aborted. 

Cause: One of the option cards trailing the CO-OCCUR control card has a misspelled word, 
or the words do not begin in the correct columns. The correct columns are 6 and 
31 for two-way tables, and 6, 31, and 56 for three-way tables. 

6. Difference between site and AFFIRMS station greater than 2,000 feet, no adjustments made. 

Cause: When using ADJUST option, elevation differences are limited to 2,000 feet. Program 
executes, but without any site adjustments. 

7. More than 5 months weather summary requested, only first 5 months processed. 

Cause: Scl£-explanatory--split data into two or three sections. 

Program RXWTHR Output Quantity 

For each station, the following output is produced: 

SUMMARY tables : One (1) page for each parameter requested. 

CO-OCCURRENCE tables : 

Two-way: One (1) page for each table per 2 months analysis; for example, a two-way 
table from May through August will have a two-page output. 

Qne (1) page for each table for each month's .analysis; for example, a three- 
way table from May through September will produce five pages of output. 
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Program RXWTHR Cost Estimates 

lor d tvpical RXl'.TilR run (one station, 5-month analysis, five summary tables, two to 
three co-occiiri'once tables) a user can expect costs and resource use as outlined below: 


Run priority 

Cost 

CAU (sec) 

Resource time (sec) 

DEMAND 

S6.00 - $7.00 

14 - 15 

22 - 25 

L 

$3.00 - $4.00 

14 - 15 

22 - 25 

M 

$2.00 - $3.00 

14 - 15 

22 - 25 

P 

$1.00 - $2.00 

14 - 15 

22 25 


RXWTHR Demonstration 


Putting together steps 1 through 4 

results in the execution 

of a very simple run of 


FIXWIDR. for demonstration purposes, a sample analysis of two stations in a single run will 
be used. Other runs are demonstrated in the Program Options section. In this example, 
Philipsburg Ranger Station and West Glacier headquai'ters are the two stations in the sample. 
Assume data file is FIREDATALIB*21-24. Remember that stations must be analyzed in ascending 
order (tiie same order in which they are retrieved from the data library). 


Step 1. r-ill out information sheets (see figs. 11 and 12). 

Step 2. Create data file for program use and execute program RXNTIIR. 

123456 7 8 

12345678901234567890123456789012345678901234567890123456789012345678901234567890 

@RUM, 


@ASG,A FIFa:DATAUB*PROGRAMS. 

OASG.A FIREDATAI,IB*21-24. 

QUSE 2., FIREDATALIB*21-24. 

@ASG,UP PHILIP. 

@USE 15., PHILIP. 

0XOT F IREDATALIB*PROGRAHS . GETDATA2 . 

24020755 24020778 

24300260 24300277 

@EOF 
@FREE 2. 

@ASG,T 2. 

@ASG,A CSSG*R1LIB. 

§XQT CSSG*R1LIB, RXWTHR 

Step 3. Transfer information from user sheet(s) to machine readable formats. 
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02 


Total \uniber Different Stations in This Run 

User’ 5 Same COLONY .5ubunit GLACIER NP 


.Unit 


NPS 


Project—'' Brief fare weather sununaries for West Glacier headquarters 


Fire leather Station Information 


Elevation 3200 

f t . , Latitude 

Fuel model- 

.Last trost — 

>ear begin 1953 

,Year end 

Site Adjustment Factors I if any): 

Aspect ( l=north, 2 = east, 3=s< 

Canopy cover 

< l=open, 2 = ' 

Duff/Soil Horizon 

Information (if 1 

Layer Duff/ Soil 

Tvoe Thickness 

1 

cm 


cm 

5 

cm 

4 

cm 

5 

cm 

Suminarv Table(s) Requested (Select 


Station name west glacier 
46 


240207 


■>/ 

.Climate class— 


.Slope class^^ 


5/ 


1978 


, Grass type - ' annual 

Date be^in 0001 .Date end 


. perenn la I 


091S 


* 

* 1_/ Use up to 80 characters 

* 

I 2/ See User's Guide, appendix C 
★ 

I y See User's Guide, appendix B 

J £/ For NFDRS indices only. See User's Guide, 

I appendix A 

J 5/ For NFDRS indices only 

* ♦ * IT* ★****♦•***★**•**★*-*★■*■ , Kkftifk'tfk-iririK'tc'k-k-kicit-Kick-M'tritk-itk-k* 


State of the weather 


X T emperature (degrees F) 

Relative humiditv 

X W ind direction (8 point) 
Wind speed (ini/h) 


_Min relative humidity i 24 h, %) 
Precip duration (last 24 h) 
^Precip amount (24 h, O.Dl in) 

_^1 hour fuel moisture 
10 hour fuel moisture Co) 


Max temperature (24 h, deg. F) 
Min temperature i24 h, deg. F) 


Max relative humidity (24 h, '’j) 


_NFDRS ERC 
_NFDRS BE 

_^Duff Moisture (24 h average, “o) 


Co-occurrence Table(s) Requested (If 2-way table desired leave last space blank. If 
selected, Wind Direction must always be listed first.) ; 


1 

WIND DIRECTION 

with 

WIND SPEED 

with 

2 

TEMPERATURE 

with 

RELATIVE HUMIDITY 

WIND SPEED 

with 

3 


with 


with 

4 


wi t h 


with 



with 


with 


Figure 11 . --RX1VTHR user information sheet for RXWTHR demonstration using the West Glacier 

weather station. 





Toial \umber Oi fferent Stations in This Run 02 


FIRE LAB 


, Unit int exp stn 


t J 

I f ““ 
Proiect— 

Demonstration of RXWTHR output 

for Philipsburg Ranger Station 


Fire leather 

Station Information 

Station 

PHILIPSBURG RS 

name 

, .Vo. 

243002 

Ele\atLon 

5280 ft., Latitude 

46 , 

21 

Climate class— 3 , 

.Slope class— 


Fuel model-- 

.Last frost 

5/ 

.Grass type — : annual_ 

. perennial 

'tear be^in 

1960 .Year end 

1977 

.Date begin 0501 

.Date end 

0931 

Site -Adjustment Factors (if any) 





Aspect 1 

[ l=north, 2=east, 5 

=south, 4=west ), Site elevation 


ft. 

Canopv cover 

( l=open. 

2=iclosed ) 





Duff/Soil Horizon Information (if Duff Moisture selected. See User^s Guide, appendix P) : 
Layer Duff/ Soil Type Thickness 
I cm 


cm 

cm 


cm 


cm 


★ ♦★**********-*****#***-»»****#*i*r*nr**T******i»n*r<»*****-<if*Hf*' 

if * 

* ★ 

J U Use up to 80 characters I 

* 

* y See User's Guide, appendix C * 

* + 

J y See User's Guide, appendix B J 

* ★ 

* £/ For NFORS indices only. See User's Guide, * 

J appendix A 

* 5/ For NFDRS Indices only 

* * ★★★♦•**i*K(f**'fr*****Hf**„it**-*#i(r*itf***'*A*'*’**'******-<*-'+** 


Summary Tablet si Requested (Select up to 15 ) : 

State of the weather 

X Temperature [degrees F) 

Relative humidit) {%) 

X w ind direction (8 point) 

Wind s])eed (mi/h) 

^iax temperature (24 h, deg. F) 

^Niin temperature (24 h, deg* F) 

relative humidity (24 h, %) 


_Min relative humidity (24 h, %) 
Precip duration (last 24 h) 
Precip amount (24 h, 0*01 in) 

_1 hour fuel moisture {%) 

10 hour fuel moisture (%) 

NFDRS ERG 
NFDRS BI 

Duff Moisture (24 h average, %) 


o~qccurrence Tdble(s) R equested (If 2-way table desired leave last s pace blank* If 
selected, Wind Direction must always be listed first*): “ 


1 WIND DIRECTION 

with 

WIND SPEED 

^with 

2 TEMPERATURE 

with 

relative humidity 

with WIND SPEED 

5 

with 


inrith 

4 

with 


with 

5 

with 


with 


Figure 12.-RX!frHR user information sheet for RXOTHR demonstration using Philipsburg weather 

Station. 


♦ J*- ♦ * n. 4- 



02 


STATION 

WEST GLACIER 240207 3200 46 3 

YEARS 1955 1978 

DATES 0601 0915 

ACTIVITY 

BRIEF FIRE WEATHER SUMMARIES FOR WEST GLACIER HEADQUARTERS 

1 2 3 4 5 6 7 

1234567890123456789012345678901234567890123456789012345678901234567890 

SUMMARY 

TEMPERATURE 

WIND DIRECTION 

END 

CO-OCCUR 

WIND DIRECTION WIND SPEED 

TEMPERATURE RELATIVE HUMIDITY WIND SPEED 

END 

RUN 

02 

STATION 

PHILIPSBURG RS 243002 5280 46 3 

YEARS 1960 1977 

DATES 0501 0931 

ACTIVITY 

DEMONSTRATION OF RXWTHR OUTPUT FOR PHILIPSBURG RANGER STATION 

RUN 

Step 4. Terminate program execution. 

QEOF 

@PIN 

RXWTHR Output Interpretation 

Identical output (format) would be produced for each station by the above demonstration, 
For brevity, only output from the second station, Philipsburg, will be discussed. 
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Sll'L'L\R'i TAB LIS 


Relative frequency of occurrence of selected for^montlily' periods 

dard deviation, median value, and value ranges are tabulated f ^ ^ j tables simply cou- 

(.ee fig. I), for parameter., wind direction r/th^moS o^ 

?ibt of relative frequency ot occurrence of each singular class t 

class. This format is displaved in figure 2. 

Interpretation of these tables should be quite straightforward; ^ 

tables of historical values. In both tables, cumulative frequencies 

addition of relative frequency values. For e.xaniple, in May (fig. 1) 85.8 percent W . p 
clt . ?6?5 perUnt , 15% percent . 12.7 percent) of the days In the sample had to.poMtnres 

Of less than 70 degrees. Similarly, only 14.2 percent of the temperature readings weic 70 
degx'ees or higher. 


When obtaining summaries of NFDRS indices the user should consider excluding tho first 21 
days of a year's record. The reason is that at least that much time is needed for t ho iucl 
moisture models to predict ^’converging" values for the 100- and 1000-hour fuels. (John E. 
Deeming, personal communicationj . 


CO-aCrURRENCL TaXBLES- -TWO-WAY 

In the two-way tables, 2 months of analysis are displayed on each page of output. One 
page (May and June) is depicted in figure 3. Values within the tables are percent frequency 
of co-occurrence of selected class values of each parameter. In May, the most frequent co- 
occurrence condition is a windspeed of 3 to 5 mi/h from the southwest (9,7 percent). Most 
\;indspeed values v;ere in the range of 3 to 8 mi/h as seen from the columnar frequency totals, 
ilus information should be used for general knowledge of an area’s weather and to screen and 
identify periods that have high probabilities of meeting potential prescription conditions,^ 
After a prescription occurrence examination by RXBURN, the information could be used in making 
new prescriptions that have high probability of occurring. 


CO-OCCURRFNCL TABLES- -TflREE-WAY 


One month of analysis is tabulated on each page of output for three-way tables. Again, 
just one page is depicted in figure 4. The last parameter entered (windspeed, in this case) 
is stratified into five value classes; the first two are stratified into 10 value classes, 

In this example, the frequency of co-occurrence of a windspeed from 6 to 11 mi/h, a temper- 
ature of 59 degrees or less, and a relative humidity of 30 to 49 percent is 9.0 percent. It 
is obtained by adding the relative frequencies of all those particular stratifications. The 
number of days in the sample is given in the bottom right portion of the table. This infor- 
mation has the same potential for use as that from the two-way tables, with tlie added flexi- 
bility of a third parameter. 


Using Program RXBURN 

To use RXBURN, the four steps designated for RXWTHR are required with a slight variation in 
the substeps. Eiie steps are relisted here with only the changes detailed. The structure of 
directive blocks is the same as for RXWTHR; they must begin with the Number of Stations card and 
end with a RUN card, 

Step 1. Fill out information sheets for RXBURN (fig. 3). 

All values up to prescription factor selections are identical to those in the RXWTHR 
sheet. A discussion on defining prescriptions for this section is given in appendix E. 

Step 2. Create data file and load program RXBURN. 

This 15 done exactly the same as in RXWTHR except the command @XQT CSSG*R1 LIB. RXWTHR is 
altered to read: 
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(axgT CSSG*R1LIB.RXBURN 

Step 3, Transfer information from user’s sheet to machine roadable formats. 

Items (a) through (d) are exactly the same as in RXWTHR, The changes in this step are 
outlined and exemplified here. In RXBURN, there is no restriction on the number of months 
analyzed in one run, 

User identification cards . --RXBURN allows the information from the second line of 
the sheet to be entered into the program for documentation. The first card of this sequence 
instructs the program that the user’s name, district, and forest follow on the next card, It 
consists of the word IDENTIFY beginning in the first column of the card. The information card 
follows with the user’s name in coliunn 1 through 15, the district’s name in 21 through 35^ and 
the forest’s name in 41 through the 55th column. 


1 2 3 4 5 6 

123456789012345678901234567890123456789012345678901234567890 

IDENTIFY 

WILLIAM SHAW TROUT LAKE RD CASCADE NP 


f. Prescription condition cards .--The first card of this sequence instructs the program how 
many prescription cards follow. The number (nn) may range from 1 to 15; the card consists of the 
word PRESCRIBE beginning in the first column, and the value nn in columns 11 and 12, right- 
justified. This card is then followed by the prescription condition cards, These cards are 
transferred directly from the user information sheet. The parameter name begins in column 6, the 
preferable RX minimum values are entered in column 33 through 35; preferable RX maximum values in 
columns 38 through 40; acceptable RX minimum values in columns 43 through 45; and acceptable maxi- 
mum values in the 48th through 50th column. All values must be right-justified. Values for wind 
direction should be entered as eight points in a clockwise direction; for example, southwest to 
northwest, northwest to northeast; and east to west. 


1 2 3 4 5 

12345678901234567890123456789012345678901234567890 


PRESCRIBE 04 

WIND SPEED 
TEMPERATURE 

DAYS SINCE LAST PRECIP 
WIND DIRECTION 


6 7 4 10 

58 59 55 63 

7 7 4 9 

NW NW W N 


g. The RUN card . --The RUN card serves the exact purpose of that in RXWTHR and must be 
the last card in each directive block. 

Step 4. Terminate program execution. 

The program is terminated the same as RXWTHR was terminated. 

In Step 3, the order in which items (b) through (f) are entered is not restricted as long 
as the order within each item is consistant with the directions stated above. This is also 
true i^^ith input for RXWTHR. 


25 


Program RXBURN Error Messages 

1. No recog, U 23 ble input cards - incorrect input stream - program aborted. 

Cause: Misspelled control word on a control card, or control card word not beginning in 
column 1. 

2. ,\t least one option card is incorrect or option word started in the wrong column - 
program aborted. 

Cause. Misspelled option afterthe PRESCRIBED control card, or option card not beginning 
in column 6. 

5. Difference between site and AFFIRMS station greater than 2,000 feet, no adjustments made. 
Cause: Self-explanatory, 

Program RXBURN Output Quantities 

for each station and/or prescription analysis in a single RXBURN run, approximately five 
pages of output are produced, with a maximum of nine in a full 12-month analysis. 

Program RXBURN Cost Estimates 

for a ti-pical RXBURN run (one station, 5 to 6 months analysis, three to four prescription 
parameters),’ a user can expect the costs and resource use outlined below. 


Run priority 

DEMAND 

Cost 

CAU (sec) 

Resource time 

$4.00 to $5.00 

8 to 9 

10 to 11 

L 

$2.00 to $2,00 

8 to 9 

10 to 11 

N 

$1.00 to $1.50 

8 to 9 

10 to 11 

P 

$0.75 to $1,00 

8 to 9 

10 to 11 


RXBURN Demonstration 

Again, combining steps 1 through 4 results in the execution of a simple run of RXBURN, In 
this demonstration, one station is used, although more could be analyzed in a manner identical to 
that shown in RXWTHR, It is again assumed that the data resides on file FIREDATALIB*21-24 , 
Remember, the weather data are read from logical unit 15, 

Step 1. Fill out RXBURN information sheet (see fig, 13). 

Step 2 . Create data file and load program RXBURN, 
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t^XBURN -- USER INF0R^4ATIQN SHEET 


Total Number Different Stations in This Run Of 

User's Name bradshaw , Subunit FIRE lab ^ Unit Int exp stn 

Projeot^^ Demonstration of RXBUPN output for Philipsburg Ranger Station 

Fire Weather Station Information - Station name philipsburg rs ^\jo, 243002 

Elevation 5280 ft,, Latitude 46 .Climate class - 3 , Slope class^;;|^ 

Fuel model^ , Last frost — ^ Grass type^'^: annual , perennial 

Year begin 1960 .Year end 1977 ,Date begin 0501 j_Date end npi 

Site Adjustment Factors fif any) : 

Aspect ( l=north, 2=^east, 3=south, 4=west ), Site elevation ft. 

Canopy cover ( l=open, 2=closed ) 

Duff/Soil Horizon Information (if Duff Moisture selected. See User's Guide, Appendix D] ; 

* * 

1 ^cm * y Use up to 80 characters * 

* 

2 ^cm * 2J See User's Guide , appendix C * 

■h * 

3 cm * V See User's Guide, appendix B * 

* 

4 cm * 4/ For NFDPS indicies only. See User's Guide, J 

+ , appendix A -v 

5 cm * V For nFDRS indices only J 

■x'k'kif^eie'k'^f Jt ★*★**+** *'W * *** •V * ★'V ♦** *«!?•>(:*★* + -Tt * * Kr + + ^ V* * * * 

Prescription Factor Selections (Check and set limits for up to IS factors.) : 

Preferable Rx Acceptable Rx 

Factor ^ti^imuln Maximum Minimum Maximum 

State of the weather 

X T emperature CUeg. F) 65 75 60 80 

X Relative humidity (%) 20 30 20 55 

^iVind direction (8 point) 

X Wind speed (mi/h) 4 9__ 0 15 

^Max temperature (24 h, deg. F) 

^Min temperature (24 h, deg, F) 

^Max relative humidity (24 h, %) 

Nlin relative humidity (24 h, %) 

^Days since last precipitation 

Precip amount (24 h, 0.01 in) 

1 hour fuel moisture (0) 

10 hour fuel moisture (%) 

^NFDRS ERG 

NFDRS BI 

^Duff moisture (2-i h average, o) ____ 

Figure 13.--RXBURN User Information Sheet for RXBURN demonstration using Philipsburg weather 

station. 



1 2 3 4 8 

1234567890123456789012345678901234567890 0 


ORUN, . . . 

@ASG,A FIREDATALIB*PROGRAMS. 

@ASG,A FIREDATALIB*21-24. 

OUSE 2. ,FIREDATALIB*21-24. 

(aASG,UP PHILIP. 

@USE 15,, PHILIP. 

SXQT FIREDATALIB*PR0GRAMS.GETDATA2 
24300360 24300277 

@EOF 
@FREE 2. 

@ASG,T 2. 

<aASG,A CSSG*R1LIB. 

@XQT CSSG*R1LIB.RXBURN 

Step 3. Transfer information from user's sheet to machine readable formats. 
01 

STATION 

PHILIPSBURG RS 243002 5280 46 3 

YEARS 1960 1977 

DATES 0510 1101 

ACTIVITY 

DEMONSTRATION OF RXBURN OUTPUT FOR PHILIPSBURG RANGER STATION 
IDENTIFY 

LARRY BRADSHAW FIRE LAB MONTANA 

PRESCRIBE 03 


temperature 

65 

75 

60 

80 

RELATIVE HUMIDITY 

20 

30 

20 

55 

WIND SPEED 

4 

9 

0 

15 


RUN 

tep 4, Terminate program execution. 
@EOP 
@FIN 
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RXBURN Output Interpretation 

Output from this demonstration is displayed in figures 5 through 8. 

SUMMARY PAGE 

Figure 5 depicts the format and values of the summary page of RXBURN output. Part 1 
consists of input infoi'mation: district name, user name, station information, sample size, and 
a prescription condition summary. Part 2 is composed of a prescription occurrence summary, Tt 
details the mean number of days per season of each prescription type, the frequency of occurrence 
for each type, and the highest 10-day period and month frequency of each prescription type. In 
this example, 8 days per year were preferable, 51 were acceptable, and 80 or 58 percent of the 
season had unacceptable conditions, July had tlic highest acceptable prescription probability at 
45 percent, and the highest 10-day period begins July 1, with a probability of 56 percent. Two 
10-day periods, June and October 1-10, had equal probabilities of preferable conditions at 10 
percent, while the least likely period for a preferable or acceptable condition to occur is the 
10-day period beginning May 1, with a probability of 22 percent of (100 percent - 78 percent). 

Figure 6 shows the seasonal progress of prescription occurrences. The sample size, mean 
number of days, total count of occurrences, and occurrence frequency are given for each prescrip- 
tion type by iO-day, monthly, and seasonal stratifications. This information should be used to 
follow the historical pattern of prescription occurrences throughout the season. It may be help- 
ful to plot the occurrence frequencies as they change over the season to obtain a better feel for 
the progression. 

Figure 7 is an example of the output from the run length summary analysis of RXBURN, The 
mean run length (number of days in a row of the same prescription type), and the values at the 
25th, 50th (median), and 75th percentile levels are given. Stratification is by 10-day and 
monthly periods. This information gives an indication of the length of time that favorable or 
unfavorable conditions tend to persist. In this demonstration, the period July 1 to July 10 
averaged a preferable run length of 2 days. At the 25th percentile level, the value was 1 day, 

2 days was the median, and the value was 3 days at the 75th percentile level. This means that 
there is a 50 percent chance of a preferable run lasting at least 2 days and a 25 percent chance 
of one lasting at least 3 days. 

The last section of output from RXBURN is the probability of a prescription tyi^e occur- 
ring in 1, 2, or 3 days, given an initial prescription type. Stratification is by month. 

Three months are displayed on each page of output. An example of this output is shown in 
figure 8, In this example, for a day in May with acceptable conditions, there is a 12 percent 
chance that tommorow^ s conditions will also be acceptable, a 9 percent chance that conditions 
2 days hence will be acceptable, and a 7 percent chance for acceptable conditions on the third 
day, If todays s conditions, in the example, are unacceptable, there is an 80 percent chance 
that they will be unacceptable the next day, and only a 2 percent chance of their being acceptable. 

System Options 

ANALYZING MORE TUAN ONE STATION 

Up to 99 stations may be analyzed in a single execution of either RXWTllR or RXBURN by 
building station data files (stations in ascending order) and "stacking'’ directive blocks (one 
for each station). This is illustrated in Step 3 of the RXWTHR demonstration. Each directive 
block must begin with a Number of Stations card and end with a RUN card, Information entered 
once within a directive block (for example, prescription condition, summary table option, or 
activity information should not be reentered in following blocks unless the information changes. 

STATIONS WITH MORE THAN 5 MONTHS WEATHER DATA 

If two or more runs of RXWTHR are required to obtain a full season analysis, set the 
Number of Stations card value to one (01) and stack two (or three) directives blocks with the 
dates of analysis cards so that 5 months are analyzed in the first run. The dates of analysis 
cards in following directive blocks should have different seasonal dates so that the entire 
season is covered in two or three runs. 
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MULTIPLE PRESCRIPTION CONDITION ANALYSIS OF ONE STATION 

Primary use of this option is to analyze two or more prescription conditions for a 
single station in one program execution. It may also be used to get more than five co- 
occurrence tables in RXWTHR. 

An initial run of the program stores computed and observed parameter values; subsequent 
runs have access to stored values, resulting m computer cost and time savings. Values are 
stored by entering a ^'T*' in Column 48 of the first directive block *s station information card. 

Stored computations are used by the entry of a ’’T" in column 46 of the second directive 
blocK^s station information cai'd (leave column 48 blank). The third and following directive 
blocks do not require' than any of the station information series be entered. There must be 
one directive block for each prescription occurrence analysis. The following example shows 
how to analyze four prescription conditions at one station. 


1 2 3 4 5 6 7 

1234567890123456789012345678901234567890123456789012345678901234567890 

01 

STATION (Steps 1 and 2 are assumed to have been completed) 

ANY STATION 123456 2345 45 H B P 3 0615 T 

YEARS 1900 1999 

DATES 0101 1215 

ACTIVITY 

DEMONSTRATION OF MULTIPLE PRESCRIPTION ANALYSIS 
IDENTIFY 

SANDRA LEE FERN R.S. ASPEN N.F. 

PRESCRIBE 01 

NFDRS ERC 2 4 1 10 

RUN 

01 

STATION 

SAME STATION 123456 2345 45 H 3 P 3 0615 T 

YEARS 1900 1999 

DATES 0101 1215 

PRESCRIBE 01 

NFDRS ERC 4 8 4 15 
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RUN 


01 

PRESCRIBE 01 

WIND SPEED 

RUN 

01 

PRESCRIBE 03 

WIND DIRECTION 
WIND SPEED 
MAX temperature 

RUN 


@EOF 

OFIN 


13 18 


S W 

4 8 

65 72 


10 20 


S NW 

2 12 

60 78 


1 2 3 4 5 6 7 

1234567890123456789012345678901234567890123456789012345678901234567890 

There is one restriction in using this option--paraineters not computed in the initial run 
cannot be used in subsequent prescription conditions, This applies only to NFDRS indices, fuel 
moisture, and duff moisture. Fuel moisture is computed when either fuel moisture or an NFDRS 
index is specified in a prescription. This also means that fuel models and site adjustments 
cannot be altered on subsequent runs. 

SITE ADJUSTMENTS 

Site adjustments factors adjust fuel moisture values (and resulting NFDRS indices) to locations 
different from those of the base weather station. This option is provided because most v^eather 
data is taken at open valley-bottom locations while much prescribed burning is done on higher 
mountain slopes. The adjustment factors are from research for Fire Behavior Officer training 
courses offered through the USDA Forest Service. Use is detailed on user information sheets and 
in table 1. A brief example of a station information card to adjust to a site 800 feet higher 
than a station, north aspect, closed canopy is shown below. Site elevations must be within 2,000 
feet of the base station. It is assumed in this example that steps 1 and 2 have been completed. 


1 2 3 4 5 6 7 

1234567890123456789012345678901234567890123456789012345678901234567890 

ANY STATION 123456 5000 45 3 3 T 1 5800 2 
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APPENDIX A - NFDRS FUEL MODELS 


t^uei MOdei A 

11 , represents western grasslands vegetated by annual grasses and forbs. 

lUbh or trees may be present, but are very sparse, occupying less than one^-third of the area. 
Lxaiupics of types where fuel model A should be used are cheatgrass and medusahead. Open 
pinyon-juiripcr, sagebrush-grass, and desert shrub associations may be appropriately assigned 
this fuel model if the woody plants meet the density criteria. The quantity and continuity 
of the ground fuels vary greatly with rainfall from year to year. 


Fuel Model B 

Mature, dense fields of brush 6 feet (1.8 m) or more in height are represented by this 
iLiol model » One-fourth or more of the aerial fuel in such stands is dead. Foliage burns 
readily. Model B fuels are potentially dangerous, fostering intense, fast -spreading fires. 
inis model is for California mixed chaparral, generally 30 years or older. The F model is 
more appropriate for pure chamise stands. The B model may also be used for the New Jersey 
pine barrens. 


Fuel Model C 

Open pine stands typify model C fuels, Perennial grasses and forbs are the primary ground 
fuel, but there is enough needle litter and branchwood present to contribute significantly 

loading. Some brush and shrubs may be present, but they are of little consequence. 
Situations covered by fuel model C are open, longleaf, slash, ponderosa, Jeffrey, and sugar 
pine stands. Some pinyon- juniper stands may qualify. 

Fuel Model D 

This fuel model Is specifically for the palmetto-gallberry understory-pine overstory 
association of the southeast coastal plains. It can also be used for the so-called "low 
pocosins" whore fuel model 0 might be too severe. This model should only be used in the south- 
east because of a high moisture of extinction characteristic. 

Fuel Model E 

Use this model after leaf fall for hardwood and mixed hardwood-conifer types where the 
hai'dwoods dominante. The fuel is primarily hai'dwood leaf litter. The oak-hickory types are 
best represented by fuel model E, but E is an acceptable choice for northern hardwoods and 
mixed forests of the southeast. In high winds, the fire danger may be underrated because 
rolling and blowing leaves are not accounted for. In the summer, after the trees have leafed 
out, fuel model E should be replaced by fuel model R. 

Fuel Model F 

Of the 1972 NFDRS fuel models, only the application of fuel model F has changed, Model 
F nov^ represents mature closed chamiso and oak brush fields of Arizona, Utah, and Colorado. It 
also applies to young, closed stands and to mature, open stands of California mixed chaparral, 
Open stands of pinyon- juniper are represented. Fire activity, however, will be overrated at 
low windspeeds and where ground fuels are sparse. 

Fuel Model G 

Fuel model G is used for dense conifer stands where there is a heavy accumulation of 
litter and domed woody material. Such stands are typically overmature and may also be suffering 
insect^ disease, wind, or ice damage- -natural events that create a very heavy buildup of dead 
material on the forest floor. The duff and litter are deep and much of the woody material 
Is more than 3 inches (7.62 cm) in diameter. The undergrowth is variable, but shrubs are 
usually restricted to openings. Types meant to be represented by fuel model G are hemlock- 
Sitka spruce, coast Douglas- fir, and windthrown or bug-killed stands of lodgepole pine and spruce. 
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Fuel Model H 


The short-needled conifers (white pines, spruces, larches, and firs) are represented by 
fuel model H, In contrast to model G fuels, fuel model H describes a healthy stand with sparse 
undergrowth and a thin layer of ground fuels. Fires in H fuels are typically slow spreading 
and are dangerous only in scattered areas where the downed woody material is concentrated. 

Fuel Model I 

Fuel model I was designed for clearcut conifer slash where the total loading of materials 
less than 6 inches (15,24 cm) in diameter exceeds 25 tons/acre (5,6 kg/m^) . After settling 
and the fines (needles and twigs) fall from the branches, fuel model I will overrate the fire 
potential. For lighter loadings of clearcut conifer slash, use fuel model J, and for light 
thinnings and partial cuts where the slash is scattered under a residual overstory, use fuel 
model K. 


Fuel Model J 

This model complements fuel model I, It is for clearcuts and heavily thinned conifer 
stands where the total loading of materials less than 6 inches (15.24 cm) in diameter is less 
than 25 tons/ acre (5,6 kg/m^). Again, as the slash ages, the fire potential will be overrated. 

Fuel Model K 

Slash fuels from light thinnings and partial cuts in conifer stands are represented by 
fuel model K. Typically the slash is scattered about under an open overstory. This model 
applies to hardwood slash and to southern pine clearcuts where the loading of all fuel is less 
than 15 tons/acre (3.36 kg/m^) . 


Fuel Model L 

This fuel model is meant to represent western grasslands vegetated by perennial grasses, 

I he principal species are coarser and the loadings heavier than those in model A fuels. Other- 
wise, the situations are very similar; shrubs and trees occupy less than one-third of the 
area. The quantity of fuel in these ai'eas is more stable from year to year. In sagebrush areas, 
fuel model T may be more appropriate. 


Fuel Model N 

This model was constructed specifically for the sawgrass prairies of south Florida. It 
may be useful in other marsh situations where the fuel is coarse and reedlike. This model 
assumes that one-third of the aerial portion of the plants is dead. Fast-spreading, intense 
fires can occur even over standing water. 


Fuel Model 0 

The 0 model applies to dense, brushlike fuels of the southeast. 0 fuels, except for a 
deep litter layer, are almost entirely living in contrast to B fuels. The foliage burns 
readily except during the active growing season. The plants are typically over 6 feet fl 8 m) 

thP Pocosins of the Virginia and 

the North and South Carolina coasts are the ideal of fuel model 0. If the plants do not meet 
the 6-foot Cl.S-m) criteria in those areas, fuel model D should be used. 


Fuel Model P 

Closed, thrifty stands of long-needled southern pines are characteristic of P fuels A 
2- to 4-inch (5- to lO-cm) layer of lightly compacted needle litter is the primary fuel ' Some 
small diameter branchwood is present, but the density of the canopy precludL more than'a 

shrubs and grass. Fuel model P has the high moisture of extinction characteristic 
of the southeast. The corresponding model for other long-needled pines is U. 
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Fuel Model Q 


Upland Alaskan black spruce is represented by fuel model Q. The stands are dense, but have 
frequent openings filled with usually inflammable shrub species. The forest floor is a deep 
layer of moss and lichens, but there is someneedle litter and small -diameter branchv\rood . The 
branches are persistent on the trees, and ground tires easily climb to the tree crows. This 
fuel model may be useful for jack pine stands in the Lake States. Ground fires are typically 
slow spreading, but a dangerous crowning potential exists. Users should be alert to such 
events and note those levels of spread component and burning index when crowing occurs. 

Fuel Model R 

This fuel model represents the hardwood areas after the canopies leaf out in the spring. 

It is provided as the off-season substitute for E. It should be used during the summer in all 
hardwood and mixed conifer-hardwood stands where more than half of the over story is deciduous. 

Fuel Model S 

Alaskan or alpine tundra on relatively well-drained sites is the S fuel. Grass and low 
shrubs are often present, but the principal fuel is a deep layer of lichens and moss. Fires 
in these fuels are 7iot fast spreading or intense, but they are difficult to extinguish. 

Fuel Model T 

The bothersome sagebrush- grass types of the Great Basin and the Intermountain West are 
characteristic of T fuels. The shrubs burn easily and are not dense enough to shade out 
grass and other herbaceous plants. The shrubs must occupy at least one-third of the site or 
the A or L fuel models should be used, Fuel model T might be used for immature scrub oak 
and desert shrub associations in the West, and the scrub oak-wire grass type in the Southeast, 

Fuel Model U 


Closed stands of western long-needled pines are covered by this model. The ground fuels 
are primarily litter and small branchwood. Grass and shrubs are precluded by the dense canopy, 
but occur in the occasional natural opening. Fuel model U should be used for ponderosa, 
Jeffrey, sugar, and red pine stands of the Lake States, Fuel model P is the corresponding 
model for southern pine plantations. 


35 



APPENDIX B - NFDRS SLOPE CLASS DEFINITIONS 


Slope class 
1 
2 

3 

4 

5 


Slope {^peroent) 
0-25 
26-40 
41-55 
56-75 


75 and above 
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NFDRS 

climate 

class 


\ 


2 


3 


4 


APPENDIX C - NFDRS CLIMATE CLASS DEFINITIONS 


1 

Thornthwaite 

humidity 

province 

Arid 


Semiarld 


Subhumid 
(rainfall 
deficient in 
summer) 


Subhumid 
(rainfall 
adequate in 
all seasons) 

Humid 


Wet 


Characteristic 

vegetation 

Desert (sparse 
grass and scat- 
tered shi'ubs) 


Steppe (short 
grass and shrubs) 


Savanna (grass- 
lands, dense 
brush and open 
conifer forests) 


Savanna (grass- 
lands and open 
hardwood forests) 


Regions 


Sonoran deserts of west Texas, 

New Mexico, southwest Nevada, 
and western Utah; and the Mojave 
Desert of California, 

The short grass prairies of the 
Great Plains; the sagebrush 
steppes and pin yon/ juniper wood- 
lands of Wyoming, Montana, Idaho, 
Colorado, Utah, Arizona, Washing- 
ton, and Oregon; and the grass 
steppes of the central valley of 
California, 

The Alaskan interior, the chap- 
arral of Colorado, Arizona, New 
Mexico, the Sierra Nevada foot- 
hills, and southern California; 
ponderosa pine woodlands of the 
West; and mountain valleys (or 
parks) of the northern and 
central Rockies, 

Bluesteni prairies and bluestem- 
oak-hickory savanna of Iowa, 
Missouri, and Illinois. 


Forests 


Rain forest 
(redwoods, and 
spruce-cedar- 
hemlock) 


Almost the entire eastern United 
States; and those higher eleva- 
tions in the West that support 
dense forests. 

Coast of northern California, 
Oregon, Washington, and southeast 
Alaska, 


^Thornthwaite 1931, 
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NFORS CLIMATE CLASS 

1 □ Arid and semiarid 

2 □ Subhumid 

3n Subhumid and humid 
4^ Wet 

Appendix C Map.— NFDRS climate class (Deeming and others 1977). 
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APPENDIX D - DUFF/SOIL HORIZON STRUCTURE 


Programs RXWTHR and RXBURN offer a newly applied duff moisture model. It must be stressed 
that the model Is theoretical and has not been validated by field tests, It is based upon 
theoretical properties of water vapor transport through porous mediums. The model was devel- 
oped by Michael Posberg [Posberg 1975) and adapted for computer usage on archival fire weather 
data by a consulting firm, it is offered strictly on an experimental basis in the present 
versions of RXWTHR and RXBURN^ Further work and field testing are planned and resulting 
changes and calibrations will be incorporated into updated versions of the programs. 

To employ this option of programs RXWTHR and RXBURN, the user needs to provide information 
about the structure of forest floor. The floor must be broken into layers of homogeneous 
materials: duff, humus, and a boundary or stable layer. At the surface, the boundary layer 
is the air, at the bottom, the layer is considered to be a sand. The user must specify the 
type of material of each layer and the thickness of each layer. The types of materials that 
the model will presently accept and the appropriate codes to be entered in the duff information 


cards are as follows: 

Duff/ soil type Code 

Lodgepole pine duff 1 

Lodgepole pine humus 2 

Ponderosa pine duff 3 

Ponderosa pine humus 4 

Douglas-fir duff 5 

Douglas-fir humus 6 

Sand 2 


The thickness of each layer is to entered by the user in centimeters. At this time the 
bottom boundary layer is assumed to be sand witli a thickness from 100 to 300 cm {approximately 
40 to 120 inches). This information still must be entered. The type entered for each layer 
is tlien used to calculate layer timelags based on physical properties of each of the layers, 
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APPENDIX E - DEFINING PRESCRIPTION CONDITIONS 

HXBtfRN offers the iiniilue feature of allowing two value ranges to be 
prescription factor in a prescription. They are defined as preferab c t - 4 . ■ ^ , 

i:rrnu,g"»;.tL»s. PrefL-.ble Lndltrons are those .deal for .ectiug f 

acceiitablo conditions are broader (and «ore often occurring) ranges of values that will still 
meet burning objectives, though not as surely as burning in preferable conditions. All values 
outside the r^ln^^e of acceptable conditions are considered unacceptable. 

As an example, a manager might plan to burn with a prescription that included a temper- 
ature of SS" or 59° F., last precipitation 7 days ago, and vnnd out of the northwest at 
6 or 7 miles per hour. Realizing this tight prescription may have a limited chance of occur- 
ring, thev consult climatological summaries from RXWTHR and decide on an acceptable prescrip- 
tion that will probably still meet burning objectives, and yet have a reasonable probability oi 
occurring. The acceptable prescription might than include a temperature of 55 to 63 h., last 

precipitation 4 to 9 days ago, and a windspeed of 4 to 10 miles per hour from the west to north, 
In terms of input to RXBURM, this prescription would look like this, 


Preferrable RX 


Acceptable RX 


Factor 


Minimum 


Maximum 


Mininiuin 


Maximum 


Temperature, 

Days since last precipitation 
Wind speed (mi/k) 


Wind direction 


Any values outside the acceptable ranges would constitute an unacceptable burning condition. 
Note that wind direction limits must be entered in a clockwise manner. 


A word of caution. Windspeed and direction can be quite unreliable except for tlie site 
at which they were collected. Consequently, the manager must rely to a large extent on 
experience and Knowledge of the effects of local topography on wind in order to interpret 
RXBURN values for these parameters. 
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APPENDIX F - RXWTHR AND RXBURN PROGRAM DOCUMENTATION 

Ihis appendix contains documentation of all routines and functions of programs RXWTHR 
and RXBURN. It is intended to serve as a reference for programers and users interested in the 
logic and structure of the pi'ograms. It provides variable identification of values stored 
in COMMON BLOCKS. 

A general overview of the structure of RXWTHR precedes its documentation. This is followed 
by a general overview of program RXBURN and documentation of the details of RXBURN that differ 
from those in RXWTHR. 

General Overview of RXWTHR 

I'rogram RXWTdlR computes summary tables and co-occurrence tables of user-selected fire 
weather parameters. All user input is read from logical unit 5 in subroutine SETUP. All 
weather data from the National Fire Weather Data Library is read in subroutine CMPUTE from logical 
unit 15. Based on user options, the pi'ogram will calculate either summary tables or co- 
occurrence tables, or both, In each run of RXWTHR, a maximum of 15 summary tables and five 
co-occurrencc tables can be computed. 

'I'he general flow of the program is as follows; 

Step 1: Read number of stations to be analyzed in mainline control prograini RXWTHR. 

Step 2; RXWTHR calls SETUP, which reads user control sequence, SETUP then calls INTERP to 
interpret user input and set program conti'ol variables and computation levels. 
INTERP then returns to SETUP, which then calculates adjustment codes if needed 
and returns to RXWTHR. 

Step 3: RXWTHR then calls CMPUTE. CMPUTE sets initial values and begins reading weather 
data and computing needed values as per user specification and summarizes data 
by calling SUM. 

Step 4: After all data have been read and values accumulated in SUM, CMPUTE calculates 
final statistics and outputs tables through two subroutines, RITEl outputs the 
summary tables, RITE2 outputs the co-occurrence tables. 

Step 5: After first station has been analyzed and queued to output, CMPUTE returns to 
RXWTHR to check for another station to analyze. If there are none the 
program terminates normally. If there is another station, general fiov^ 
starts at Step 2 again until all stations have been analyzed. 

In the following pages, the COMMON BLOCKS used by RXWTflR are first defined. Mainline 
program RXW'illR is then documented. All other routines and functions are then documented in 
alptmbetical order, 

RXWTHR Common Blocks 

Program RXWI’IIR uses nine labeled COMMON BLOCKS and no unlabeled COMMON BLOCKS, These are 
defined as follows: 

CO^^MON /FIRE/ : Contains parameters used in system fire model routine. 

VALUES SET : Values in /FIRE/ are set in subroutine FUELS and subroutine CURING, all 

variables are REAL. 
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TRIABLE 


D 

OD 

OOD 

GOOD 

DP 

PTH 

WOOD 

ERB 

GHERB 

moD 

lOOD 

31D 

:iOD 

DOOD 

3RB 

HRBC 

)FC 

)OD 

} 

:RnC 

WON /LABELS/ ; 


DESCRIPTION 


Weight of dead and l-hour fuels (pounds per square foot) 

Weight of dead and down 10-hour fuels (pounds per square foot) 
Weight of dead and do™ 100-hour fuels (pounds per square foot) 
Weight of dead and down Ij 000-hour fuels (pounds per square foot) 
Weight of dead and dovm 1-hour fuels + cured green fuels 
Indicates first call to subroutine ERCSC 
Depth of fuel bed (inches) 

Height of dead fuels (inches) 

Height of live fuels (inches) 

Moisture content of woody fuels (percent) 

Moisture content of herbs (percent) 

Surface area to volume ratio lierbaceous fuels 

Surface area to volume ratio woody fuels 

Surface area to volume ratio 100-hour dead and down 

Surface area to volume ratio 1-hour dead and dom 

Surface area to volume ratio 10-hour dead and down 

Surface area to volume ratio 1,000-hour dead and down 

Weight of herbaceous material 

Weight of cured herbacious material 

Wind factor 

Weight of woody fuels 

Moisture of extinction of dead and down fuels (percent) 

Moisture content of cured herbaceous material (percent) 


Contains alphanumeric labels for output. All variables are integer and 
set in BLOCK DATA, 
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AY : 

nnscRiPTioN 



.n't, 17) : 

Header labels for all parameters 



.2(9) : 

Header labels for wind direction 



.3(10) : 

Table labels for summary tables 



vK17) : 

Units labels for all parameters 



L5(17) : 

Short header labels for co-occurrence tables 



vlTl[(12) : 

Month labels (JAN^-DEC) 



^IMON /LIM/ 

Contains class level limits for all parameters. All values in /LIM/ are 
integer and are set in BLOCK DATA except for LOW, which is set in the out- 
put routines. 

ray : DESCRIPTION 



MIT (9, 17) 

Contains class value limits for all parameters except STATE OF THE WEATHER 

MSIV(IO) 

Contains class value limits for STATE OF THE WEATHER 



)WC9) 

Contains lower' class limits foi’ output tables 



)NLMON /ODSERV/ 

Contains values from weather observation, computed moisture 
other values needed for computation 

valuesj 

and 

IRAYS 

: DESCRIPTION 



DilC22) 

: Contains weather observation, plus computed values 



OBCl) 

; Station number 



ODC2) 

: Year 



OIiC3) 

: Month 



OB(4) 

: Day 



OB(5) 

; State of the weather 



OB (6) 

: Temperature 
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I0B(7] 

I0B(8) 

rOB(9) 

IOB(103 

IOB(ii:) 

I0B(12) 

I0B(14) 

rOB(15) 

I0BC16) 

I0BC17) 

I0B(18) 

I0B(19) 

IOB(20) 

I0BC21) 

IOB(22) 

STRLfCC5,9) 

STRUCCI,!) 
STRUCCI,2) 
STRUCCI, 3) 
STRUCCI, 4) 
STRUG (I, 5) 
STRUCCI, 6) 
STRUCCI, 7) 
STRUCCI, 8) 
STRUCCI, 9) 

lELV 
I DUFF 


Relative humidity 
Duff moisture (computed) 

Wind direction 
Windspeed 

Stick moisture (observed) 

Maximum temperature 
Minimum temperature 
Maximum relative humidity 
Minimum relative humidity 
Precipitation duration 
Precipitation amount 
Moisture index (1, 2, or 3) 

1-hour fuel moisture (computed) 

10-hour fuel moisture (computed or observed) 

1978 NFDRS ERC (computed) 

1978 NFDRS BI (computed) 

Contains duff/soil horizon structure information as input from SETUP, 
are up to five horizons. 

Soil or duff type code 

Thickness of layer (centimeters) 

Particle temperature timelag (seconds) 

Particle moisture timelag (seconds) 

Bulk density layer (grams per cubic centimeter) 

Particle density (grams per cubic centimeter) 

Hydraulic conductivity (centimeters per second) 

Moisture coefficient (dimensionless) 

Temperature coefficient (dimensionless) 

Adjustment elevation code (1--3) 

Duff moisture option specification (0 - no, 1 = yes) 


There 
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IDAYB 

IDAYE 

ISLO 

lYEARB 

lYEARE 

ELEV 

LAX 

MONTHB 
MONTI IE 
REPEAT 

SAVE 

NUMSTA 

MKT 

MNT 

MXRII 

MNRH 

I APT 

CO^t^!ON /STATS/ 


ARRAY 

MEAN(5,3,153 

SDC5,3,15] 

IXTREMC5, 3, 15) 
XKYEAR(5,3,15) 


Day to begin sample inclusion 
Day to end sample inclusion 
Adjustment slope code (1 or 2) 

Year to begin sample 
Year to end sample 
Station elevation 
Station latitude (12) 

Month to begin sample 
Month to end sample 

Logical (LI) variable that specifies if run is to read from disk file 
created by SAVE option 

Logical (LI) variable that specifies if disk file is to be created for 
subsequent REPEAT runs 

6-digit station identifier code 

Maximum temperature - ALPHA format 

Minimum temperature - ALPHA format 

Maximum relative humidity - ALPHA format 

Minimum relative humidity - ALPHA format 

Precipitation amount - ALPHA format 


Contains statistical summary values for summary and co-occurrence tables, 

Arrays for use in summary tables are stratified by month (5); 10-day 
period of month (3); and summary table number (15). Arrays for co-occurrenco 
are stratified by table number (5); month (5); first order class number (iO^ ; 
second order class number (10); and third order class number (5) transferr'ccl 
into a vector. 


DESCRIPTION 

Accumulates values stratified as described above for summary tables. 

Accumulates sums of squares and then computed values of standard deviations 
are entered, stratified as described above. 

Holds maximum values of parameters for use in range display in summary tables 
Holds minimum values of parameters for use in range display in summary tables 
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MEDIAN(2250) : Holds values ordered by class for computation of median values for display 
of summary tables. Values are stratified by period of month (3), month (5)^ 
summary table number (IS), and class interval number (12), tranformed into a 
vector. 

CUMCO (12500) : Cumulator array for co-occurrence tables, stratified as described above. 

COUNT (5, 3) : Holds sample size stratified by month and period of month. 

Program RXWTHR 

Narrative: Program RXWTHR is the mainline control program for the computation of clima- 
tological summary co-occurrence tables. It keeps track of the number of stations to be analyzed 
and the number that have been analyzed. It resets the program or terminates it according to the 
number of stations ret^uested, 

The logic is; 

Step 1: Set run number to zero. 

Step 2: Read number of stations to be analyzed. 

Step 3: If at EOF, terminate run normally; if not, continue to Step 4. 

Step 4: Call SETUP 
Step 5; Call CMPUTE 

Step 6: Rewind disk file created by SAVE option (if used) and add 1 to the run number. 

Step 7; Check to see if last station has been analyzed. If not, return to Step 2, If so, 
rewind weather data file, set run number to zero again, and return to Step 2 to see 
if stations are to be reanalyzed with perhaps different fuel models, slope classes, 
or co-occurrence tables. 

COMMON BLOCKS USED: NONE 

A listing of all SUBROUTINES and FUNCTIONS that compose the body of RXWTHR now follows. 

They are listed in alphabetical order and follow the same format as the mainline program docu- 
mentation of RXWTHR. 

RXWTHR Ructions and Subroutines 

SUBROUTINE ADSORP(I) 

Narrative: Subroutine ADSORP is used in duff moisture calculations during periods that 
experienced rainfall. ADSORP is called by PERCOL, which was called by DUFF, which was called 
CMPUTE. The parameter I is passed into ADSORP and is the current layer in the duff horizon 
bhat is being simulated. 


lOMMON BLOCKS used: /BL0CK2/, /OBSERV/ 


""he logic is: 


■tep 1: Solve for effective percolation rate PPP, and then set QQQ and RRR from it. 
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Step 2; Call ROOT- -ROOT returns value for TilETA after solving a third degree polynomial equation. 
Step 3: Compute particle moisture content adsorption rate (percent moisture content per second). 
Step 4: Compute actual moisture content change (amount = rate x time). 

Step 5; Adjust THETA (volumetric content of voids) by amount adsorbed by particles, 

Step 6: Calculate effective hydraulic conductivity due to water adsorption. 

Step 7: Return to PERCOL. 

SUBROUTINE ALPHA(RF,AF) 

Narrative: Subroutine ALPHA is passed the array RF, an array containing real values. It 

ENCODES the values and if the value is zero, sets it equal to an alpha blank in an A4 format. 

This makes all output of relative frequency values alpha, for easier readability. AF is the 
array that is passed back for output, ALPllA is called by RITEl and uses no COMMON BLOCKS. 

BLOCK DATA 

Narrative: Routine BLOCK DATA sets data constants for all LABELED COMMON BLOCKS that are 

given values through the use of DATA statements. 


COMMON BLOCKS USED: /OPSET/, /LABELS/, /LIM/ 


FUNCTION CALV(F) 

Narrative: This real function converts temperature in degrees F to temperature in degrees 
Kelvin and returns to calling program. CALV is a function from the NFDRS, 

SUBROUTINE CMPUTE 

Narrative; This routine is the major computation section of the program. It reads the 
weather data and calls all other subroutines in the process, It computes statistical summaries 
and directs output procedures by calling output routines at appropriate places in the flow. 

COMMON BLOCKS USED: /SET/, /OBSORV/, /STATS/, /OPSET/, /LABELS/, /FIRE/, /LIM/ 


The logic is: 

Step 1: Initialize all constants, including site adjustment factors (lADJUS), zero out all 
values in COMMON BLOCK /STATS/. 

Step 2; If fuel model is needed as per computation level (L) , match the fuel model that was 
requested (MOD) with the test array of (MODI) and get fuel bed parameters by calling 
FUEL(II). Then set the day of last frost (LFROST) to a month and day value for use 
in subroutine CURING. 

Step 3: Compute fire weather station pressure, adding 200 feet to the station elevation if 
Alaskan. 

Step 4: Begin DO LOOP of reading weather data from logical unit 15, 


47 



step 5: 


Step 6: 


Step 

7: 

Step 

8: 

Step 

9: 

Step 

10: 

Step 

11: 

Step 

12: 

Step 

13: 

Step 

14: 

Step 

15: 

Step 

16: 

Step 

17: 

Step 

18: 

Step 

19: 

Step 

20: 

Step 

21: 

Step 

22: 


If REPEAT is run, skip calculations and read values from disk file. Otherwise ^ make 
checks for correct station and to see if the day just read is within the time period 
requested. If maximum/ minimum values are requested for analysis, but missing from 
data record, read another record (Step 4) . 

Determine if a new season has started. A new season is said to have started if the year 
changes or there is a 30-day break in the records. If new season has started, set the 
logical variable INIT to .TRUE. This allows for initial values to be set. 

Call processing option setting routine OPTFIX. 

Call data testing routine TESTWX. If there is an error in the data, write a message 
and read the next record. 

If DUFF MOISTURE is requested, CALL DUFF (INIT, IFIRST) 

If no fuel moistures or NFDRS indices are requested, go to Step 20. 

Compute solar declination by calling function DECL and passing month and day. 

Get 1- and 10-hour fuel moistures from subroutine ONETEN. 

Make site adjustments if requested and add to values of 1- and 10-hour fuel moistures. 
If no NFDRS indices are requested go to Step 19. 

Compute 100- and 1,000-hour fuel moistures by use of routines MlOO or MIOOA, depending 
on value returned from OPTFIX. 

Call subroutine CURING to determine state of green vegetation. 

Call subroutine ERCSC to determine the day's NFDRS indices. 

Assign computed values of ERC and BI to elements 21 and 22 of lOB, respectively. 

Assign correct values of 1- and 10-hour fuel moisture to elements 19 and 20 of lOB, 
respectively. 

T£ computations are to be saved for next run use, write them to the disk file that is 
called 'tape 2'. 

Call subroutine SUM, This routine does the intermediate statistical summaries.’ 

Read a new record (Step 4). 


At this point all records for the station have been read and summed in subroutine SUM and are ready 
for final computations and output. 


Step 30; If no summary tables have been requested, skip to co-occurrence table output section. 
Otherwise, proceed. 

Step 31; Begin DO LOOP for each summary table requested (1 to NMEANS) . 


Step 52: Get summary parameter type, order, and write out header by call to RITEl . If fuel 

model was used, write out fuel model letter and slope class. Then write out table- tops. 

.Step 33; Begin DO LOOP for each month of data (1 to NMONTH) . 
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step 34: 
Step 35: 

Step 36 ; 
Step 37: 

Step 38: 

Step 39: 
Step 40: 
Step 41: 

Step 42: 


Begin DO LOOP for each period of the month (1 to 3), 

Set cumulative frequencies, relative frequencies, extreme values, and if cumulative 
frequency has reached 50’’^ get median value (MED) by calling subroutine MIDPT2 (this is 
for first class interval). 

Begin DO LOOP for each class interval beginnin at 2. (Number 1 was set in Step 35), 

Set cumulative frequencies, etc,, as in 35 till median value is determined. Then end 
Step 36 LOOP. 

Compute period mean, standard deviation, and output period mean, standard deviation, 
median, range values, and relative frequency distribution by calling RITEl. 

Continue LOOP from Step 34. 

Continue LOOP from Step 33, 

Repeat process, this time by monthly stratifications by summing three periods for each 
month, and outputting by call to RITEl. 

Repeat entire j^rocess for each summary table requested (Step 31) 


At this point all summary tables have been computed and written to the p'^int queue. Flow continues 
in computation and output of co-occurrence tables. 


Stop 

50: 

Step 

51: 

step 

52: 

Stop 

53: 

Stop 

54: 

Stop 

55: 

step 

56: 

Step 

57: 

Step 

58: 

Step 

59: 

Step 

50: 

Step 

61: 

Step 

62; 

Step 

63: 


If no co-occurrence tables were requested, return to RXWTHR. 

Begin DO LOOP for number of co-occurrence tables requested. 

Set control points based upon options and order of options. 

If this table is a 2-way table, skip to Step 70. Otherwise, proceed. 

Begin DO LOOP for number of months 

Sum total number of days in all periods of the month and write out header, fuel model, 
and slope class if needed, 

Zero out row and column totals and begin computation loops, 

First write out table-tops by calling subroutine RITE2. 

Begin DO LOOP for number of rows (class intervals) for first parameter in the table 
(1 to L4 [Option 1,1]). 

Compute an entire row of percent frequency of co-occurrence for the class values and 
output the entire row by calling RITE2, 

Continue Step 58 LOOP. 

Write out column frequency totals by calling RITE2 (there are two sets). 

Compute total frequency (row and column totals should be 100%) and write totals out 
along with the total count for the month (number of days in the sample). 

Repeat for each month in data set (Step 54) , 
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step 64: Repeat foi' each table (Step 51). 


The next block is for computation and output of 2-way tables of co-occurrence. 


Step 70 
Step 71 
Step 72 
Step lit 
Step 74 
Step 75; 
Step 76 
Step 77; 


Begin DO LOOP for number of months. 

Write headers, fuel model, and slope class if needed. 

Zero out row and column arrays (ROWTOT and COLTOT) . 

Write out table-tops. 

Compute simple 2-way co-occurrence frequencies and output results one row at a time. 
Compute row and column total, sample site, and write it out. 

Repeat for each month (Step 70). 

Repeat for each table (Step 51). 


At this point all summary and co-occurrcnce tables have been computed and queued to the printer. 


Step 80: If computations were saved on disk, put end of file mark (hOh) at the end of the data. 
Step 81: Return to RXWTllR, 

SUBROUTINE CURING(INIT) 

Narrative: Subi-outine CURING Is designed to model tJie pi'ogrossLon of green vegetation 
through the growing season from dead, through gveenup, througli transition, and again to tlie 
dead or cured state. The subroutine was developed by the NlMlRS work unit, and taken in entirety 
from the NFDRS. Full documentation of this subroutine Is available through appropriate NFDRS 
publications. 


COMMON BLOCKS USED: /OBSERV/, /FIRE/ 


FUNCTION DECL(MONTii,IDAY) 

Narrative: This function computes the solar declination angle and is taken from NFDRS 
research. It is used for calculations of lOO-hour fuel moisture. 


FUNCTION DLIGHTCLAT,DECL) 

Narrative: This function computes the total hours of dayliglit for tlio station latitude 
and the solar declination angle as computed in DliCL. It Is also used In calculations of 100- 
hour fuel moisture. DLIGIIT was derived from NFDRS algorithms. ' 

SUBROUTINE DUFPCINIT.IFIRST) 

Narrative: This subroutine controls the calculation of duff moisture. It Is still in the 
developmental stages and should be used with caution, "Tljc model was developed from numerical and 
analytical solutions of the diffusion forms of the mass continuity equation and the first law of 
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thermodynamics. Analytical solutions provided a fundamental framework to evaluate nonlinear 
interactions obtained in the numerical solutions. Dimensional analysis was used to define the 
relationship between the soil properties used in the model,'’ (Fosberg 1975.) 


COMMON BLOCKS USED; /OBSERV/, /BLOCKl/, /BL0CK2/ 


Tile logic is: 


Step 1: Set subroutine constants. 


Step 2: If it is the first time the subroutine has been called on this run, set up physical 
characteristics of the duff horizon. These factors include (by layer): 

1. Packing ratio (BETA) 

2. Porosity (PHI) 

3. Vapor Diffusivity (NU) 

4. Thermal Diffusivity (KAPPA) 

5. Layer Moisture Timelag (TAURHO) 

6. Layer Temperature Timelag (TAUTMP) 

7. Modified Exponential Exchange Functions (temperature [CIIITMP]) 

B, Modified Exponential Exchange Functions (vapor [CHIRHO]) 

9. Depth of layer (DEPTH). 

Step 3: [f it is initial season, call PROFLl to set up initial temperature and relative humidity 
profiles of the horizons. 


Step 4: Begin simulation. 

Step 5: Get input values of temperature and relative humidity. 

Step 6: Begin 24-hour duff moisture simulation (2-hcur increments), 

Stop 7: Use harmonic analysis to model temperature and relative humidity cycles for each time 
increment. 

Step 8: If TPTAMT is greater than zero, and simulation time is less than precipitation duration 
it is still raining. GO TO Step 10. 


Step 9: This means that it has either stopped raining or did not ram. If ENG is less than 

PIBSAT, regular solution will be used and SWITCH is .TRUE, for the layer, Otherwise, 

moisture is lost at the above saturation point rate as defined by DTHETA, This is done 
for each layer, and skips to Step 11. 

stop 10- Call PERCOL. This sets a percolation rate of free water through the horizons which is 
^ adsorped to an extent by a value of DTHETA, and contributes to the overall moisture 

content of the horizon by a volumetric amount that is determined by solution of a third 
SS^rorpolyLi.! .quatiL by subroutine ROOT. This yaluo is ,iote™inod for each layer, 
and what is adsorped by one layer is not available for adsorption by the next. 

Step 11: Determine vapor pressure for each layer and set top layer boundary conditions and inter 

mediate layers to last value of layer above. 


Step 12: Modify layer temp and density by exchange functions. 

Step 13: Recalculate vapor pressure using new layer temperature (sat v.p.). 


Step 14: Compute current saturation density. 
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Step 15: If EMC greater than FIBSAT, current density - saturation density. 

Step 16: Compute vapor pressure and layer relative humidity from equation of state. 

Step 17: If Sh^ITCEl of layer is .TRUE, get particle moisture by calling DUFFMC which calculates 

EQMC of particles based on temperature and relative humidity. Compute average^ simulate 
next time period, until 24 hours have been simulated, then return to CMPUTE with I0B(8) 
duff moisture, 

SUBROUTINE DUF FMC (TEMP, HUMID, EQM, TYPE) 

Narrative: Subroutine DUFFMC computes duff particle equilibrium moisture contents based 
on a regression shown in the Wood Products Handbook, using predictors of temperature and 
relative humidity. It also computes EQMC^s for soil or loam boundary layers. 

FUNCTION EQMC (TEMP, RH) 

Narrative: This function provides the same information as DUFFMC, but is in the form of 
a function that does not have the capability for soil and boundary layer computations, 

SUROUTINE ERCSC(ERC, BI) 

Narrative: This computational routine determines potential fire behavior from the NFDRS 
fuel model and the current fire weather. It then computes the fire danger indices based on 
the potential fire behavior. The model is based on work of Rothermel (1972) and Albini (1976) 
and modified by Deeming and others (1977), This routine is taken directly from the NFDRS 
system and again is fully documented there. 


COMMON BLOCKS USED: /FIRE/ 


SUBROUTINE ERROR(I) 

Narrative: This subroutine is called by several routines in the system and is passed 
the value I, The routine is called when an eri-or is detected. An error message is printed 
by means of a computed GO TO statement and processing continues or the program is aborted, 
depending on the severity of the error. 

FUNCTION FREEZE ( ITMN , NUNfDAY) 

Narrative: This function determines if a killing frost has occurred and is called by 
CURING, A killing frost is considered to have occurred if the minimum temperature has been 
less than 33 degrees for a total of 5 days, or if a minimum temperature of 25 degrees or less 
occurs, This logical function is taken from the NFDRS. 

SUBROUTINE FUEL (II) 

Narrative; This routine sets the parameters of the fuel bed for analysis or modeling by 
ERCSG, The routine contains the 20 NFDRS stylized fuel models, and assigns the correct values 
to the variables in COMMON BLOCK /FIRE/. FUEL is called by CMPUTE, 

I BROUT I NE 1 NTERP (M VA LUE , MOPT I N) 

Narrative: This routine interprets what was read in SETUP, sets program directives, and 
omputation levels (L) . It is passed the alphanumeric words that were read by SETUP, matches 
-hem with its word dictionary, and sets option parameters and computation levels. If no 
latch is found, ERROR is called with the appropraite value passed. 


OMMON BLOCKS USED: /SET/, /OPSET/ 
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FUNCTION lOlJT(KSET) 


Narrative: This function examines the option number and sets a value (TOUT) for use in 
COMPUTED GO TO in the output routines. 

FUNCTION IRND(X} 

Narrative: This integer function rounds a real number then performs a FORTRAN standard 
IF IX operation on it, 

INTEGER FUNCTION ISUBl (I , J, K, L) 

Narrative; This function converts a four-way stratification into a vector for computer 
storage on a computer that is restricted to three dimensional arrays. This particular function 
IS used in computing median values in the summary array MEDIAN. 

INTEGER FUNCTION ISUB2 (I , J, K, L,M) 

Narrative: This function converts a five-way stratification into a vector for computer 
storage. This particular function is used in co-occurrence stratification of array CUMCO. 

INTEGER FUNCTION JET (I) 

Narrative; This function groups precipitation and ERG values into categories of graduated 
value, 

SUBROUTINE MlOO (DECL, LAT, BNDRY) 

Narrative: The routine is called by CMPUTE to calculate 100-hour fuel moistures, It is 
the preferred method of 100-hour calculations. Maximum and minimum temperatures and relative 
humidities are used to calculate extreme values of equilibrium moisture contents which are then 
integrated with DLIGIIT resulting in a weighted average value that is used in the 100-hour 
computations. This routine was developed by the NFDRS work unit and includes BNDRY values. 


COMMON BLOCKS USED: /OBSERV/ 


SUBROUTINE Ml OOA (IWTFLG, LAT, DECL, BNDRY) 

Narrative: This routine developed by the NFDRS work unit is for computing values of 100- 
hour fuel moisture when maximum and minimum temperatures and relative humidity are not avail- 
able, It uses default values of moisture recovery at night based on climate classes and 
yesterday’s 100-hour fuel moisture content. It then computes a boundary condition that is 
used in computations of 1,000-hour fuel moistures. 


COMMON BLOCKS USED: /OBSERV/ 

SUBROUTINE MIOOO (BNDRY, MClOOO, INIT) 

Narrative: This NFDRS routine computes 1,000-hour fuel moisture by using a 7-day average 
boundary (BNDRY) condition as a driver function in the calculations. The value for BNDRY is 
computed in either MlOO or MIOOA. 

SUBROUTINE MIDPT2 (M,C, X, CO, J, L, LI) 

Narrative: This subroutine computes the median value of an order ranked distribution by 
methods for large sample sizes described by Nie and others (1975). M is the computed median 
value, C is the cumulative frequency, X is the last class's relative frequency, CO is the 
sample size, J is an on/off switch, and L and LI are upper and lower class value limits. 
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SUBROUTINE MODE CRF , J, ISET) 

Narrative: This subroutine selects the modal class value from a relative frequency 
distribution, RF is the distribution* J is an on/off switch, and ISET is an option delimiter 
set in FUNCTION lOUT. 

INTEGER FUNCTION NDEX (I , J* K, L,M) 

Narrative: This function computes the vector element of a 5-way stratification adjust- 
ment vector. The stratifications are (1) aspect, (2) elevation, (3) slope, C^) canopy cover, 
and (5) month, 

SUBROUTINE ONETEN 

Narrative: This function computes 1- and 10-hour fuel moisture values from the daily 
fire weather observation. This routine was developed by the NFDRS and has several additional 
data integrity checks for use in this program. 

The logic is: 

Step 1: Set temperature and relative humidity adjustment factors to compute fuel bed tempera-* 
ture and relative humidity based on cloud cover. 

Step 2; Calculate equilibrium moisture content by FUNCTION EQMC, 

Step 3; If stick moisture is reported and in allowable limits, use it and calculate 1-hour 
fuel moisture from 10- hour value and return. 

Step 4: Use regression equation to compute 1- and 10-hour values if stick moisture is not 
reported. 

Step 5: Return to CMPUTE. 


COMMON BLOCKS USED: /OBSERV/ 


SUBROUTINE OPTFTX 

Narrative: This NFDRS routine computes which processing option will be used for fuel 
moisture computations. lOPTIN == 1 is the preferred option and lOPTIN is the least desirable 
option (due to mimimal data), 


COMMON BLOCKS USED: /OBSERV/ 

SUBROUTINE PERCOL 

Narrative: This subrouitne is called by DUFF, in the case of precipitation, to compute 
the percolation rate of free water through the duff horizons. 

COMMON BLOCKS USED: /OBSERV/, /BL0CK2/ 


The logic is: 

Step 1; Compute precipitation rate. 
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Step 2: Set percolation rate as the minimum value of the precipitation rate of the hydraulic 
conductivity of the layer. 

Step 3: Call ADSORP to compute THETA and DTHETA for the first layer. 

Step 4: Adjust EMC and DMC as affected by THETA and DELTAM set SWITCH to FALSE. 

Step 5: Repeat for lower layers. 

Step 6: Return to DUFF. 

SUBROUTINE PROFLl 

Narrative: This subroutine is called by DUFF to initialize soil/duff horizons to initial 
season values. It is called at the beginning of each new season. The bottom layer (organic) 
and the lowest horizon (nonorganic^ are assumed to be at fiber saturation. A linear gradient 
is established from the bottom layer to the surface relative humidity value. 

CO^^MON BLOCKS USED: /OBSERV/, /BL0CK2/ 

The logic is: 

Step 1: Set boundary layer initial conditions. 

Step 2: Determine bottom layer subscript. 

Step 3: Assume bottom layer (organic) and lowest horizon (nonorganic) are saturated. 

Step 4: Establish linear gradients of temperature and humidity between top and bottom layers. 
Step 5: Calculate vapor pressure and vapor densities of each horizon. 

Step 6: Return to DUFF, 

SUBROUTINE RELHUM 

Narrative: This adaptation of an NFDRS routine computes the relative humidity given wet 
and dry bulb temperatures. 

SUBROUTINE RITEl 

Narrative: This is the main output subroutine for summary tables. It first calls ALPHA 
to change all output to a ”A’' format and then goes to the correct writing sequence via a series 
of computed GO TO statements operating on variables passed in the calling routine CMPUTE, 

COMMON BLOCKS USED: /LABELS/, /SET/, /LIM/ 

SUBROUTINE RITE2 

Narrative: This is the main output subroutine for co-occurrence tables. It first calls 
ALPHA to convert REAL values to an ''alpha*’ format. It then goes to the proper entry point via 
computed GO TO statements operating on variables passed in the calling routine CMPUTE, 

COMMON BLOCKS USED: /LABELS/, /LIM/ 
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SUBROUTINE R00T(NC, III, JJJ, KKK, PPP, QQQ, RRR,ADFM) 

Narrative: This routine solves an Nth degree polynomial by numeric methods. The value 
returned in ADFM(l) is the volumetric moisture content the duff void space for solution to 
free water due to rain infilitration through the layers, The routine was given to the project 
by Michael Fosberg, 

SUBROUTINE SETUP 

Narrative: The call to SETUP reads user run information cards and ’’sets up'’ the run. It 
is called by RXIVTHR, It calls INTERP and returns to RXWTHR. By reading an interpreting control 
card word, the program executes proper reading sequences by means of computed GO TO statements. 

COMMON BLOCKS USED: /OPSET/, /SET/, /OBSERV/ 

The logic is: 

Step 1: Read control card word. 

Step 2: Interpret word from internal word file. If RUN card, call INTERP, otherwise execute 
proper read sequence and return to step 1. 

Step 3: Compute adjustment indices if required. 

Step 4: Return to RXWTHR, 

SUBROUTINE SUM 

Narrative; This is an intermediate computational routine for accumulating parameters 
needed for means, standard deviation, range, median, and frequency distribution calculations. 

COMMON BLOCKS USED: /SET/, /STATS/, /OBSERV/, /OPSET/ 

The logic is: 

Step 1: Compute month index and period number and increase value in sample size array COUNT 
by 1, 

Step 2: Begin DO LOOP for each summary table requested. 

Step 3: Set array pointer value K, based on summary parameter. 

Step 4: Accumulate values for mean and sum of square computations. 

Step 5: Compute interval class of summary parameter. 

Step 6: Increment correct vector element of MEDIAN by 1. 

Step 7: Set high and low values for range analysis, 

Step 8; Continue summary DO LOOP until executed once for each summary parameter. 

Step 9: If any co-occux'rence tables requested, begin DO LOOP from 1 to NTABLE. Otherwise 
return to CMPUTE, 

Step 10; Set class interval and index values for each variable in current table. 
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Step 11: Increment correct element of co-occurrence vector CUMCO by 1. 

Step 12: Continue DO LOOP. 

Step 13: Return to CMPUTE, 

SUBROin^INE TESTWX 

Narrative: This NFDRS routine tests each day's weather data for errors, computes relative 
humidity of need, and checks or computes precipitation duration. It sets some potentially 
undefined variables in the lOB array when used on REPEAT runs, 

COMMON BLOCKS USED: /OBSERV/ 


FUNCTION U20FNC(SGBRT,U20,WNDFC) 

Narrative; This function reduces the 20-foot observed windspeed to wind at wind flame height 
for calculations in ERCSC, It is called from ERCSC and is fuel bed dependent, 

SUBROUT INE VAPOR (T , EWS) 

Narrative: This subroutine computes saturation vapor pressure using formulas from the 
Smithsonian tables (Smithsonian Institution 1939). It is from the NFDRS and returns pressure 
in millibars. 

General Overview of RXBURN 

Program RXBURN computes local prescription condition frequencies based on user-defined 
prescription conditions. All user input (except for the number of stations to be analyzed) 
is read in SETUP from logical unit 5. All weather data are read from logical unit 15 in sub- 
routine CMPUTE, All output (except for error messages) is done from subroutine MEAN. RXBURN 
allows up to 15 prescription factors to be entered into a prescription. 

The general flow of the program is as follows: 

Step 1; Read number of station to be analyzed in mainline control program RXBURN. 

Step 2: RXBURN calls SETUP, which reads user control sequence, SETUP then calls subroutine 
INTERP to interpret user input, set program control variables and computation levels. 
INTERP then returns to SETUP, which then computes adjustment codes if needed and returns 
to RXBURN. 

Step 3: RXBURN then calls CMPUTE. CMPUTE sets initial values and begins reading weather data 
and computing needed values as per user specification. 

Step 4: CMPUTE then calls SUMMER, which tests to see if the prescription factors just read/ 

computed constitute a preferable condiiton, an acceptable condition, or an unacceptable 
condition. SUMMER then accumulates values needed for mean values, run lengths, and 
persistence probabilities. SU^^MER then returns to CMPUTE, which reads the next day's 
weather , 

Step 5: When all weather data has been read, analyzed, and summed, CMPUTE calls MEAN, Sub- 
routine MEAN computes all final statistics and writes all output to the printer queue. 
MEAN then returns to CMPUTE. 

Step 6: CMPUTE returns to RXBURN to check for another station to analyze, If there is none, 
the program terminates normally. If there is another station, general flow returns 
to step 2. This is repeated until all stations have been analyzed. 
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In the following pages, the COMMON BLOCKS used by RXBURN are first defined. Mainline 
program RXBURN is then documented, followed by all other subroutines and functions in alpha- 
betical order. 

RXBURN Common Blocks 

Program RXBURN uses seven labeled COMMON BLOCKS and no unlabeled COMMON BLOCKS. Ml values 
are set in the programs and none in DATA statements. The COMMON BLOCKS are defined as follows: 

COMMON /BLOCKl/ ; This common block is the same as /BLOCKl/ of RXiVTIIR. 

COMMON /BLOCK2/ : This common block is the same as /BLOCK2/ of RXWTIIR. 

COMMON /FIRE/ : ^This common block is the same as /FIRE/ of RXWTIIR. 

COMMON /OBSERV/ : This common block is the same as /OBSERV/ of RXWTIIR, 

COMMON /SET/ : Values needed for the setting up and control of Program RXBURN are held 

in this common block. It is basically a combination of values from 
/OPSET/ and /SET/ of RXWTHR, plus variables needed in RXBURN. 

VALUES SET ; Most values in this block are set in subroutine SETUP and subroutine 

INTERP. 


ARRAYS : DESCRIPTION 

ACT(20] : Contains alphanumeric activity information (20A4 elements). 

DIST(4) : Contains alphanumeric name of forest district (4A4) . 

FOREST(4) : Contains alphanumeric name of forest (4A4) . 

IRX0P(15) : Contains pointer values for extracting values from the array lOB, 

LAMAT(15, 2, 2) : Contains alphanumeric prescription limits stratified by prescription param- 

eter, preferable prescription (RX) minimum limit, preferable RX maximum 
limit, acceptable RX minimum limit, and acceptable RX maximum limit, 

LIMIT (15, 2, 2) : Contain integer limits for prescription, stratified as in LAMAT, 

NAME(43 : Contains alphanumeric user name (4A4) , 

NAMSTA(4) : Contains alphanumeric station name (4A4) . 

NRX0P(15) : Contains the order values stored in IRXOP. 


VARIABLES 

ADJUST 

ELEV 

lASP 

ICOV 

IDAYB 


DESCRIPTION 

Logical switch for inclusion of adjustment option, 
Elevation of base fire weather station (feet) , 

Aspect adjustment code for use when ADJUST is true. 

Canopy cover adjustment code for use when ADJUST is true. 
Day of year to begin inclusion of days in sample, 
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IDAYE 


Day of year to end inclusion of days in sample. 

IDUFF : On/off switch that indicates if duff moisture is to be calculated. 

Program RXBURN 

Narrative: Program RXBURN is the mainline control program for the assessment of local 
prescription condition frequencies. It keeps track of the number of stations to be analyzed 
and the number that have been analyzed. 

Tlie logic is: 

Step 1; Set run number to zero. 

Step 2; Read number of stations to be analyzed, 

Step 3: If at EOF, terminate the run normally. If not, continue to Step 4. 

Step 4: Call SETUP. 

Step 5: Call CMPUTE, 

Step 6: Rewind disk file created by SAVE option (if used) and add 1 to the run number. 

Step 7; Check to see if last station has been analyzed. If not, return to Step 2. If it has 

been analyzed, rewind weather data file, set run number to zero, and return to Step 2. 

SUBROUTINE ADSORP 

Narrative: Same as RXWTIIR. 

FUNCTION CALV 

Narrative: Same as RXWTHR. 

SUBROUTINE CMPUTE 

Narrative: This routine is a major weather reader and computational block. It performs 
the exact processes that make up the balance of CMPUTE in RXWTHR, except that no statistics 
are computed. In RXBURN, the statistics and output are handled in subroutines SUMMER and MEAN. 

COMMON BLOCKS USED: /SET/, /OBSERV/, /FIRE/ 

The logic is: 

Steps 1 through 22 are essentially the same as RXWTHR except for two items. First, alpha input 
LIMIT values are DECODED to integer format, and secondly, SUMMER is called instead of SUM at 
Step 21. 

Step 30: Call SUMMER to close out last observation strings, and ENCODE LIMIT back to alpha 
formats for output. 

Step 31; Call MEAN and put EOF on logical unit 2 (disk file). 

Step 32; Return to RXBURN. 
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RXBURN Functions and Subroutines 


SUBROUTINE CURING 

Narrative: Same as RXWTIIR* 

FUNCTION DECL 

Narrative: Same as RXIVTHR, 

FUNCTION DLIGIfr 

Narrative: Same as RXWTHR, 

SUBROUTINE DUFF 

Narrative: Same as RXWTHR. 

SUBROUTINE DUFFMC 

Narrative; Same as RXWTHR. 

FUNCTION EQMC 

Narrative; Same as RXWTHR. 

SUBROUTINE ERROR 

Narrative: Same as RXWTHR. 

FUNCTION FREEZE 

Narrative: Same as RXWTHR. 

SUBROUTINE FUEL 

Narrative: Same as RXWTHR. 

FUNCTION INDEX 

Narrative; Same as RXWTHR. 

SUBROUTINE INTERP 

Narrative: Same as RXWTHR with changes for RXBURN, 
FUNCTION IRND 

Narrative: Same as RXIVTHR. 

SUBROUTINE Ml 00 

Narrative: Same as RXWTHR, 

SUBROUTINE MIOOA 

Narrative: Same as RXWTHR. 

SUBROUTINE Ml 000 

Narrative: Same as RXWTHR. 
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SUBROUTINE MAX 


Narrative: This subroutine selects the maximum occurrence frequency for summary output 
of highest month and 10 days of frequency of prescription occurrence. 

SUBROUTINE MID 

Narrative: This subroutine selects quartile values of run length as called by MEAN at 
the right time* and sets a switch that indicates value has been computed. 

SUBROUTINE MEAN 

Narrative: This subroutine makes the final calculations of accumulation arrays and does 
all program output, 


COMMON BLOCKS USED; /SET/, /OBSERV/, /STATS/, /STATS2/ 

The logic is: 

Step 1: Set initial values and call ZER02, 

Step 2: Set total counts for seasonal, monthly, and 10-day period values, and write out first 
page head and prescription entries. 

Step 3: Compute average number of days of each prescription ty^e per season and compute 
associated frequencies and output values. 

Step 4: Select highest monthly frequency of each prescription type and write them out, As many 
as two highest months are allowed for cases of equal frequencies. 

Step 5: Select 10-day period oP highest prescription frequency of each type and output, 

Again, up to two periods are allowed. 

Step 6; Write out monthly and 10-day summaries of seasonal progression of prescription 
condition occurrence. 

Step 7: Write out seasonal totals. 

Step 8: Calculate and output mean values and quartile values of run length by monthly and 
iO-day periods. 

Stop 9: Using a simple Markov model, compute and write persistence and transition probabil- 
ities. Stratification is by month. 

Step 10: Return to CMPUTE. 

FUNCTION MTU 

Narrative: This function sets a month pointer for use in output of the alpha month. 

SUBROUTINE ONETEN 

Narrative: Same as RXWTHR. 

SUBROUTINE OPTFIX 

Narrative: Same as RXWTHR. 
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SUB[10UTINE PERCOL 


Narrative: Same as RXlvTUR. 

SUBROUTINE PROELl 

Narrative; Same as RXWTHR. 

subroutine RELUM 

Narrative: Same as RXWTHR. 

SUBROUTINE ROOT 

Narrative: Same as RXWTHR, 

SUBROUTINE SET 

Narrative: Sets switch values for use in determination of maximum frequencies in MEAN. 
SUBROUTINE SETUP 

Narrative: Same as SETUP in RXOTHR with appropriate changes for RXBURN values. 

SUBROUTINE SUMMER 

Narrative: This routine fills the accumulating arrays contained in COMMON BLOCK /STATS/. 

The logic is: 

Step 1: Set initial values--month number, period number, and increment sample size counter. 

Step 2: Check for beginning of season (INIT) and set appropriate value as needed. 

Step 3: Determine if day is "preferrable," "acceptable," or "unacceptable." 

Step 4: Increment mean value accumulator arrays. 

Step 5: Increment appropriate run length counters. 

Step 6: Determine if prescription "type" string was broken. If not go to Step 8. 

Step 7: Set appropriate run length counter arrays depending on what "type" of string was broken. 
Step 8: Set counting arrays for use in Markov model. 

Step 9. Set all today's values to holding value to become yesterday's value on the next day, 

Step 10: Return to CMPUTE, 

SUBROUTINE -TESTWX 

Narrative: Same as RXin’HR, 

FUNCTION U20FNC 

Narrative: Same as RXWTHR. 
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SUBROUTINE VAPOR 


Narrative: Same as RXWTHR. 


SUBROUTINE ZCROI 

Narrative: Sets all values in COMMON BLOCK /STATS/ to zero for each run of RXBURN, 


COMMON BLOCKS USED: /STATS/ 

SUBROUTINE ZER02 

Narrative: Sets ail values in COMMON BLOCK /STATS2/ to zero for each run of RXBURN, 
COMMON BLOCKS USED: /STATS2/ 
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The Iiueriiiountain Station, heaclquaiteird ui O/'deii, 
Utah, is one of eight legional e.xpciiuieiu stations ciuiig.ed 
with piovidmg seicntific knowledge to help icsoueoe 
manageus meet human needs anil protect foicst and langc 
ecosystems. 

The Intermountain Station includes the States of 
Montana, Idaho, Utah, Nevada, and western Wyoining. 
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aie maintained in: 


Boise, Idaho 
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State University) 
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University of Montana) /,'■ L; 
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with Biiaham Young 


Provo, Utah (in cooper 
University) 

Reno, Nevada (in cooperation with the Univcisity 
of Nevada) 




